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Abstract. The disadvantage of 'traditional' multidimensional continued fraction al- 
gorithms is that it is not known whether they provide simultaneous rational approxi- 
mations for generic vectors. Following ideas of Dani, Lagarias and Kleinbock-Margulis 
we describe a simple algorithm based on the dynamics of flows on the homogeneous 
space SL(d, Z)\ SL(d, R) (the space of lattices of covolume one) that indeed yields best 
possible approximations to any irrational vector. The algorithm is ideally suited for a 
number of dynamical applications that involve small divisor problems. We explicitely 
construct renormalization schemes for (a) the linearization of vector fields on tori of 
arbitrary dimension and (b) the construction of invariant tori for Hamiltonian systems. 
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I. Introduction 

This paper has two main goals. The first one is to introduce a new multidimensional 
continued fraction algorithm that is ideally suited for different dynamical applications. 
The algorithm can be used in order to effectively deal with small divisors whenever 
quasi-periodicity with several frequencies is an essential feature of a problem. Our 
second goal is to demonstrate the strength of the algorithm by applying it to KAM 
theory. We use a renormalization approach to prove two theorems of KAM type. The 
method, being conceptually very simple, is also very general, and allows us to consider a 
wide class of frequency vectors. For reasons of clarity we restrict our attention to vectors 
satisfying an explicit Diophantine condition (valid for a set of vectors of full Lebesgue 
measure); generalizations to more general frequency vectors follow straightforwardly 
from the approach presented here, and will be detailed in a separate publication. 

The classical continued fraction algorithm produces, for every irrational a G M., a 
sequence of rational numbers p n /q n that approximate a up to an error of order \jq\. 
The first objective of this paper is to develop a multidimensional analogue that allows 
us to approximate any irrational ex. G by rational vectors. The drawback of all 

'traditional' multidimensional continued fraction algorithms is that it is not even known 
whether they provide simultaneous rational approximations for almost Lebesgue almost 
all ex. In the case d > 4 the only result in this direction is a recent computer- assisted 
proof of the almost everywhere strong convergence for ordered Jacobi- Perron algorithm 
^OlE]- However, even in this case an explicit description of the set of bad vectors seems 
difficult. For example, the existence of 'noble' vectors, that is vectors corresponding to 
a periodic continued-fraction expansion, for which approximations do not converge, is 
rather unsatisfactory. 

The algorithm we employ here does not suffer from such pathologies. Following 
Lagarias' seminal ideas in [T3] . our approach is based on the dynamics of the geodesic 1 
flow on the homogeneous space T\G with G = SL(d, M.) and T = SL(d, Z). Notice that 
F\G may be identified with the space of lattices in M. d of covolume one or, equivalently, 
with the Teichmuller space of flat (^-dimensional tori. 

The problem of multidimensional continued fractions may be formulated in the fol- 
lowing way. Given a vector ex G R d_1 find a sequence of matrices G GL(d,Z), 
n G N, such that the 'cocycle' corresponding to the products = T'Mj'C™- 1 ) . . . j^ 1 ) 
exponentially contracts in the direction of the vector u) = ( f ) G M. d and exponentially 
expands in all other directions. Thus, the cocycle should have one negative Lyapunov 
exponent and d—1 positive. In this spirit, our algorithm comprises the following steps: 

(1) With every ol G R^ 1 associate the orbit {C(t) : t > 0} C SL(d,R), where 

<w - (V ?) { ld -f A) • <"> 

1 The term 'geodesic' is slightly inaccurate when d > 2. The orbits of the flow on T\G we are 
discussing here correspond in fact to geodesies on the unit cotangent bundle of the space T\G/ SO(g0 
for a certain family of initial conditions. Only for d — 2 the cotangent bundle can be identified with 

r\G. 
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ld-i denotes the (d — 1) x (d — 1) unit matrix. 

(2) Fix a fundamental domain T of SL(d, Z) in SL(d, R). Given a sequence of 
times ti < t 2 < ■ ■ ■ — > oo, use classical reduction theory to find the matrices 
P (n) G SL(d, Z) that map the points C{t n ) to T. 

(3) Define the nth continued fraction map by 



a (-D _> a w = iii^ ±3* n .21 



where a^ ) = a: and 



rp(n) _ I J- U L 



(n) + (n)' 
12 

T.(n) ,(n) 
6 21 fc 22 



GSL(d,Z) (1.3) 



is the nth transfer matrix defined by = y( n )p(™ x \ 

Since the action of SL(d, R) on R d_1 by fractional linear transformation defines a group 
action, we have 

P (n) rv -J- n (n) 

qW = " +P "„ (1-4) 

T P 2 1 « +P22 

where 

p(n) = T {n) T {n~l) . . . ^(1) = Pl2 g Z ). ( L 5) 

V P21 P22 ) 

Dani [5 and Kleinbock-Margulis ^2] observed that Diophantine properties of a trans- 
late to divergence properties of the corresponding orbit {TC(t) : t > 0} C T\G in the 
cusps of F\G. We exploit these results to show that, under mild Diophantine condi- 
tions on a. (satisfied by a set of a. of full Lebesgue measure, cf. Sec. 12. 6J) . there is a 
sequence of times £„ such that the transfer matrices are uniformly hyperbolic in a 
sense made precise in Sec. 12.81 This fact allows us to develop renormalization schemes 
for vector fields and Hamiltonian flows that had previously been constructed only in 
dimension one [21] or for very special choices of a El 120] • I n particular, we 

obtain renormalization-based proofs of the following theorems. 

We say us G M. d is Diophantine if there are constants e > 0, C > such that 

||fc||( d - 1 )( 1 + e )|fc. a ,| > c, (1.6) 

for all k G Z d - {0}. Note that we may assume without loss of generality that oj is 
of the form uj = ( ™) with a G R d_1 . Condition (jl.6|) then translates to a standard 
Diophantine condition on a, see Sec. (12 .6|) for details. 

Theorem 1.1. For any real analytic vector field v on T d , d>2, sufficiently close to a 
constant vector field with Diophantine vector u) G R d ; there is e > 0, an analytic curve 
(— e, e) 3 s i— > p s G R d and an analytic conjugacy h isotopic to the identity between the 
flow generated by v + p s and the linear flow 4>t(x) = x + t(l + s)cc» on T d , t > 0, for 
each \s\ < e. Moreover, the maps v i— > h and v i— > p are analytic. 

Let us emphasize that the result holds for all analytic vector fields close to a constant 
one without any additional conditions, such as preservation of volume etc. The second 
theorem deals with the case of Hamiltonian vector fields. Let B C R d , d > 2, be an 
open set containing the origin, and let H° be a real-analytic Hamiltonian function 

H°(x, y)=u-y+ l - T yQy, (x, y)eT d x B, (1.7) 
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with uj G M. d and a real symmetric d x d matrix Q. It is said to be non-degenerate if 
det Q ^ 0. 

Theorem 1.2. Suppose H° is non-degenerate and uj is Diophantine. If H is a real 
analytic Hamiltonian on T d x B sufficiently close to H°, then the Hamiltonian flow of 
H leaves invariant a Lagrangian d- dimensional torus where it is analytically conjugated 
to the linear flow 4>t{ x ) — x + tuj on T d , t > 0. The conjugacy depends analytically on 
H. 

Outline of the proof of Theorem \l.l\ Consider a vector field X(x, y) = uj + f(x) + y 
where x G T d , and y G M. d is an auxiliary parameter. The vector field f(x) is a 
sufficiently small analytic perturbation of a constant vector field. We may furthermore 
assume that uj = ( ™) for some Diophantine a G M> d ~ l ; this achieved by a rescaling 
of time. The aim is to find a value of parameter y = y^ such that the vector field 
X{x, y w ) is linearizable to a constant vector field identically equal to uj by means of an 
analytic transformation of the coordinates on T d . 

Renormalization is an iterative process, and we thus assume that after the (n — l)th 
renormalization step the vector field is of the slightly more general form 

X n ^(x, y) = oj {n -^ + / n _!(x, y) (1.8) 

where a/* -1 ) = ( ) and 

a (n i) j g gj ven by continued fraction algorithm, cf. 
()1.4p . The Fourier modes of f n -\ are smaller than in the previous step, and decay 
exponentially as ||fe|| — > +oo. We define a cone of resonant modes by a relation = 
{k£Z d : \k- uj( n ~V\ < a n _ 1 ||fc||}. The nth step requires the following operations. 

(1) Eliminate all Fourier modes outside of the resonant cone I^-i- 

(2) Apply a linear operator corresponding to a coordinate transformation given by 
the inverse transfer matrix 

(3) Rescale time to ensure that the frequency vector is of the form oj^ = ( a< f ) ). 

The conjugate action on the Fourier modes is given by k i— > T T^ 1 k. It follows from 
the hyperbolicity of that this transformation contracts for k G if cr n _i is small 
enough. This gives a significant improvement of the analyticity domain which results 
in the decrease of the estimates for the corresponding Fourier modes. As a result, 
all Fourier modes apart from the zero modes get smaller. To decrease the size of the 
latter, we choose a parameter y = y n in such a way that the corresponding zero modes 
vanish, and then consider a neighbourhood of y- values centred at y n . That is, the 
auxiliary parameter y is used to eliminate an instability in the direction of constant 
vector fields. To get enough control on the parameter dependence we perform an affine 
rescaling of this parameter on every renormalization step. One can then show that the 
corresponding sequence of parameter domains is nested and converges to a single point 
y = y w for which the initial vector field is indeed linearizable. 

In order for the scheme to be effective, the sequence of stopping times t n and the 
sizes of the resonant cones defined by the sequence of a n must be chosen properly. 
Large intervals St n = t n — t n _i improve hyperbolicity but, on the other hand, worsen 
estimates for the norms ||T^ n ^||, \\T^ \\. Similarly, if a n is too small the elimination 
of non-resonant modes will give large contributions; on the other hand, for large values 
of a n the multiplication by T T^ n+1 ) will not yield a contraction for k G I+. As we 
shall show below a right choice of sequences {(t n , a n )} can be made, depending on the 
Diophantine properties of the vector uj. □ 
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The idea of renormalization was introduced to the theory of dynamical systems by 
Feigenbaum 6J in the late 1970's. In the case of Hamiltonian systems with two degrees 
of freedom MacKay proposed in the early 1980 's a renormalization scheme for the con- 
struction of KAM invariant tori (221 • The scheme was realized for the construction of 
invariant curves for two-dimensional conservative maps of the cylinder. An important 
feature of MacKay' s approach is the analysis of both smooth KAM invariant curves 
and so-called critical curves corresponding to critical values of a parameter above which 
invariant curves no longer exist. From the point of view of renormalization theory the 
KAM curves correspond to a trivial linear fixed point for the renormalization transfor- 
mations, while critical curves give rise to very complicated fixed points with nontrivial 
critical behavior. MacKay's renormalization scheme was carried out only for a small 
class of Diophantine rotation numbers with periodic continued fraction expansion (such 
as the golden mean). Khanin and Sinai studied a different renormalization scheme for 
general Diophantine rotation numbers |14j . Both of the above early approaches were 
based on renormalization for maps or their generating functions. Essentially, the renor- 
malization transformations are defined in the space of pairs of mappings which, being 
iterates of the same map, commute with each other. These commutativity conditions 
cause difficult technical problems, and led MacKay [23] to propose the development 
of alternative renormalization schemes acting directly on vector fields. The same idea 
was realized by Koch who proves a KAM type result for analytic perturbations 
of linear Hamiltonians H°(x,y) = uj ■ y, for frequencies ui which are eigenvectors of 
hyperbolic matrices in SL(d, Z) with only one unstable direction. Notice that the set of 
such frequencies has zero Lebesgue measure and in the case d = 2 corresponds to vec- 
tors with a quadratic irrational slope. Further improvements and applications of Koch's 
techniques appeared in [U H3 EOl Ell E| , emphasizing the connection between KAM and 
renormalization theories. The results of this paper illustrate that such a programme 
can indeed be carried out in considerable generality. Another direction was followed in 
[18J, presenting a computer-assisted proof of the existence of MacKay's golden mean 
critical renormalization fixed point in the context of Hamiltonian vector fields with two 
degrees of freedom. 

Other renormalization ideas have appeared in the context of the stability of invariant 
tori for nearly integrable Hamiltonian systems (see e.g. [HHUE] and references therein). 
Inspired by quantum field theory and an analogy with KAM theory, this approach uses 
an iterative resummation of Poincare's Lindstedt series to prove its convergence. 

For the sake of transparency we have restricted our attention to Diophantine vectors 
uj. A more detailed analysis under weaker Diophantine conditions is, in principle, 
possible within the present framework. It is however a fundamental open problem to 
state a sharp (i.e. the weakest possible) Diophantine condition under which the above 
conjugacy can be established. The answer to this question is known only in the classical 
case d = 2 j2Hj where the Diophantine condition is of Brjuno type. 

It would also be interesting to see whether the multidimensional continued fraction 
algorithm presented here will allow generalizations of other one-dimensional renormal- 
ization constructions. A concrete challenge is for instance the extension of the recent 
results on the reducibility of cocycles over irrational rotations by Avila and Krikorian 


In the next section we introduce the multidimensional continued fraction algorithm, 
and include a discussion of its hyperbolicity properties required in the renormalization 
schemes. The remainder of the paper provides a detailed account of two exemplary 
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cases, the renormalization of vector fields (Section EJ) and Hamiltonian flows (Section 

ED- 

2. Multidimensional continued fractions and flows on homogeneous 

SPACES 

2.1. Flows on homogeneous spaces. Let us set G = SL(d, 1R) and T = SL(<i, Z), and 
define the diagonal subgroup {E l : t G E} in G, where 

E t = diag(e ri V . . , e rd ') (2.1) 

with constants satisfying the conditions 

d 

ri,...,r d „! < < r d , = 0. (2.2) 

i=i 

The right action of E l on the homogeneous space T\GT\G generates the flow 

$* : Y\G -> TM ^ TME 1 . (2.3) 

Since G is a simple Lie group with finite center, is ergodic and mixing [2*4*j . 

Let JF C G be a fundamental domain of the left action of T on G. Recall that, by 
definition of the fundamental domain of a free group action, 

(J PT = G, T n P.F = for all P £ T - {1}, (2.4) 
Per 

and hence, for any given M £ G, there is a unique family of P(t) e T such that 

M(t) := P{t)ME l G ^ (2.5) 

holds for all t G E. 

2.2. A convenient parametrization. Let us consider those M G G which can be 
written as 

'1 a\ ( A 0" 



M =(o ljU 7 J (2-6) 

where A G Matd_i(M) (the space of real id — 1) x id — 1) matrices), a,/3 G ]R d_1 are 
column vectors, 7 G K with 7 > 0. This yields a local parametrization of G for the set 

G+ := {(0 l) \% 7) G G : A G Mat rf-i( R )> Me K d_1 , 7 e »>o} , (2.7) 

which is particularly convenient for our purposes. All other matrices are either of the 
above form with 7 < instead, or may be written as 

M = S <*) ( * °) (2.8) 



JD I J \ '/3 7^ 

where S G V is a suitably chosen "signed permutation matrix", i.e., every row and every 
column contains one and only one non-zero coefficient, which is either 1 or —1. In the 
following we will stay clear of the parameter singularity at 7 = 0, and thus may assume 
without loss of generality S — 1. 

To work out the action of a general element T G G in the above parametrization, 
consider 

I : M — M : TM (2.9) 

where 

r (210) 
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M is as above and 



J) I) \ T (3 7 
A short calculation yields the fractional linear action 

Tucx + t l2 

a i-> cx = -= , (2-12) 

and 

7 >-> 7 = ( T *2i« + *22)7, (2-13) 
and more complicated expressions for A, /3 which will not be needed in the following. 

2.3. Multidimensional continued fractions. Let t Q = < t\ < t 2 < . . . —>■ oo be 

sequence of times, with gaps 

6t n :=t n -t n -! (2.14) 

chosen large enough so that P(t n ) ^ P(i n _i), where P(t) is defined by (12. 5|) . The se- 
quence P*" - -* := P(t n ) of matrices in T may be viewed as the continued fraction approxi- 
mants of the vector ct, which are the best possible for suitable choices of a fundamental 
domain T and times t n , see jT^j. Let us furthermore put := M(t n ) with M(t) as 
in (|2.5j) . and define by the decomposition (|2.fij) . i.e., by 

m< "'=(j a r)(v») 7°>)' 

From iW» = P^MP> and (f2~T2j) . (gUSj) we deduce 

aW = " +P " (2-16) 

T P21 « +P22 

and 

7 M = ( T p£Wpg)e^ 7 (2.17) 

where 

/ p(") „W\ 

^ = T ( n ) ^ ■ (2-18) 

V P21 P22/ 

It is evident that if the components of ( T a, 1) are linearly independent over Q, then 
7 7^ implies 7 ( - n * ) 7^ for all n > 0. 

We shall later employ the transfer matrices defined by P^ = TwpC" -1 ). Here, 

M (n) = T (n) M (n-l) E St n implies 

T (n) (n-1) , 

aW = ^ (2.19) 

and 

7 (») = (^a^- 1 ) + 4^)e^ 7 ("- 1 ), (2.20) 

where 

TW = I t5») h ] • (2 - 21] 

t 21 ^22 
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2.4. Siegel sets. In dimensions d > 2 it is difficult to describe the geometry of a fun- 
damental domain T . To overcome this problem, C. Siegel introduced simply connected 
sets Sd C G which have the property that they contain T and are contained in a finite 
number of translates PJ 7 , P G T. Consider the Iwasawa decomposition 



M = nak 



where 



n 



n u 



\ 



12 



and k G SO(cf), with Uij, cij G M, 
Sd — {nak : n G J 7 ^, aj > 



Uld \ 

Ud-X,d 

1 J 

dj > 0, CL\ 

>/3 



/ai 



V 

1. Then 



\ 



>0(j 



,d - 1), k G SO(d)} 



(2.22) 



(2.23) 



(2.24) 



is an example of a Siegel set [25J; here jF/y denotes a compact fundamental region of 
(T n N)\N, where N is the upper triangular group of elements of the form n as above. 



2.5. Dani's correspondence. We assume from now on that ri, . . . , r^-i = —1, = 
d — 1, i.e., 

E t = diag(e-*, . . . , e~*, e (d - 1}t ). (2.25) 

Let us denote by | ■ | the maximum norm in IR d_1 . A vector a G M^ 1 is called 6<26% 
approximable or of bounded type, if one of the following equivalent conditions is satisfied. 

(i) There exists a constant C > such that 

\ka + mf-^kl > C (2.26) 

for all m G Z d-1 , fc G Z - {0}. 

(ii) There exists a constant C > swc/i t/jai 

Iml^lm- a + &| > C (2.27) 

/or a// m G Z^ 1 - {0}, jfe G Z. 

The statements (i) and (ii) are equivalent in view of Khintchine's transference prin- 
ciple ([1] Chapter V). 

We recall Dani's correspondence in the following proposition (cf. [5], Theorem 2.20). 

Proposition 2.1. JTie orte {TME t : i > 0} ; with M as in (f2~Kf). zs bounded in Y\G if 
and only if the vector ct is of bounded type. 

The reason why the parameters A, (3, 7 are irrelevant in the statement is that the 
family of matrices 

W(t)=E- t ^ (3 (2.28) 

is bounded in G for alH > 0. 

The boundedness of the orbit {TME t : t > 0} implies of course that there is a 
compact set C G G such that M(t) G C for all t > 0, with M(t) as in (|23J) . 
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2.6. Diophantine conditions. A vector a G M^ 1 is called Diophantine, if there exist 
constants e > 0, C > such that 

\m\ {d - 1)il+t) \m -a + k\>C (2.29) 

for all m G Z d_1 — {0}, k G Z. It is well known that Diophantine vectors form a set of 
full Lebesgue measure @]. 

Let us show that (|2.29J) implies the inequality 



|| fe ||(d-l)(l+ e )| fc . w j >c , (2_ 30 ) 

for all fe G Z d - {0}, where w = (?), cf. (0)1 . With k = {f), $FM) yields 

|fe • w| > Clm]-^- 1 )^ > Clfel-^^ > CUfcll-^- 1 ^ 1 ^ (2.31) 

for all m G Z d_1 — {0}, k G Z. In the case when m = 0, we have k ^ (since A; 7^ 0) 
and thus ()2.3()j) holds trivially. 2 Note also that ()2.3()j) evidently implies ( )2.29j) . however 
with different choices for C in both inequalities. 

Following jTH] we define the following function on G, 

5(M)= inf | T fcM|. (2.32) 
fcez d -{o} 

It is easily checked that 5(M) is invariant under left action of T, and may thus be viewed 
as a function on T\G. In terms of the Iwasawa parametrization (|2.22j) and the Siegel 
set Sd defined in ()2.24|) we have the following estimate. 

Lemma 2.2. For M = nak G S d as in ()2.22j) . (|2.24|) . there are constants < C\ < C 2 
such that for all < a d < 1 

da d < 8{M) < C 2 a d . (2.33) 

Proof. Since ||je|| <C |x| <C ||a;|| for all x G M. d , we may prove the statement of the 
lemma for the function 

SAM) = inf || T fcM|| (2.34) 
fcez d -{o} 

instead. 3 Due to the rotational invariance of the Euclidean distance we may assume 
that k G S0(d) is the identity. 

Proof by induction. The statement trivially holds for d — 1. Therefore let us assume 
the assertion is true for dimension d — 1. The jth coefficient of the vector T kM is 

( T kM )j = (kj + a r (2.35) 

^ i=i ' 

Since a\ — > 00 when — * 0, this implies that when taking the infimum in ()2.32j) we 
must take k\ — for all sufficiently small a d . Thus we now need to estimate 



inf max 

fc e Z<i-i-{0}2<j<d 



I kj + ki Ua ) a, = a-, 1 d inf max f kj + L Uj,- ) a,- 

\ / fe eZ <i-l- { 0}2< J -<d V ^ V 



(2.36) 



2 Note that every admissible constant in l|2.29[l needs to satisfy C < 1/2; to see this, choose m = 
(1,0, ...,0), and fc € Z such that |«i + fc| < 1/2. 

3 In the following, A -C B means 'there is a constant C > such that A < CB\ If A « B « i wc 
will also use the notation Ax B. 
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where T k = (k 2 , . . . ,kd), % = a\^ d a,j so that a 2 ■ ■ • a d — 1. Now 

i-i 



inf max 

fc e Z<i-i-{0} 2 <J< d 



( % + fci Ui i ) 1 



where M = na with 



/l «23 • • ■ U 2 d \ 



n 



\ 



Ud-l.d 

1 / 



: 5 d -i(M) ^ 5 d -i(M) 
(a 2 \ 



(2.37) 



V 



(2.38) 



It is easily checked that M e <S d _i, so by the induction hypothesis, for suitable constants 
< Cw_i < C 2 ,d_i, we have 



Ci,d_i a d < 5d_i(M) < C 2 ,d-i a d , 



(2.39) 



provided = < 1. So for a d sufficiently small and a l /^a d < 1, we have 

Ci jd _i a d < 5 d (M) < C 2id _! o d . (2.40) 

In the remaining case a d > 1, all aj are bounded from above and below by positive 
constants, and hence 5 d _i(M) is bounded from above and below by positive constants. 
Furthermore a d > 1 implies a^ 1 ^ < a d , and, in view of our choice of the Siegel set, 



a-2 ■ ■ ■ a-d 3> a d So 



a d « a- 1 ^ < 



(2.41) 
□ 



and the required bound follows from ()2.3(ij) also for the case a d > 1. 

Lemma 2.3. Choose M as in ()2.6|) and suppose cx satisfies condition ()2.29|) . TTien 
i/jere ernste a constant C > smc/i i/iot /or allt >0 

5(ME t ) > C'e- 9t (2.42) 
(d-l)e 



# = 

d+(d-l)e' 

Proof. Let us put T fc = ( m. k) with m e and fceZ. Then 

( T me-*, ( T ma + A;)e (d - 1)t ) W(t) \ 



[2A3) 



5(ME l ) = inf 

(m,fc)€Z d -{0} 



> inf |( T me-',( T ma + )t)e^ 1) ')| (2.44) 

(m,fc)eZ d -{0} 

since W(i), as defined in (|2.28p . is bounded in G for all t > 0. Furthermore for t 
sufficiently large 

inf \( T me- t ,{ T ma + k)e {d - 1)t )\= inf |( W,( T ma + )t)e (d - 1)t )| 

(m,fc)eZ d -{0} meZd-i—fO}, fceZ 

(2.45) 

which, in view of the Diophantine condition ([2.29)1 . is bounded from below by 



> 



inf 



T me-* Clml-^^^e^^ 



-{0} 



e~ et inf K^Clajl-^^^ 

meZ d - 1 -{0} 



where x = e^ 9 ^'m. We conclude the proof by noting that 



(2.46) 



inf |( T aj,q a j|-^- 1 ^ 1+e ))| > inf |( T y, C\y\^ d - 1){1+e) )\ > 0. (2.47) 



meZ d - 1 -{0} 



yGK d - 1 -{0} 
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□ 

The fact that e = implies 9 = is consistent with Dani's correspondence. On the 
other hand, 9 < 1 for any e < oo. 



2.7. Norm estimates. Let 



| denote the usual matrix norm 
||M|| :=bu P MM. 

£C^0 \\ X \\ 



(2.48) 



Proposition 2.4. Choose M = M^°> as in ()2.6jl . and suppose ex. satisfies condition 
()2.29|) . Then there are constants c\, c 2 , C3, C4, C5, C6 > stzc/i that for all n e N U {0} 

||M (n) || < ciexp[(d-l)W n ], (2.49) 

< c 2 exp(0* n ), (2.50) 

< c 3 exp[(d0 + l-0)t n ], (2.51) 

< c 4 exp[(d-l + 0)t n ], (2.52) 

< c s exp[(l-9)5t n + d9t n ], (2.53) 

< c 6 exp[(d-l){l-9)8t n + d9t n ). (2.54) 
Proof. For any M £ 5j as in ()2.22|) we have, for all < < 1, 

||M|| <ax = (a 2 ---a d )" 1 < a; ( ^ 1} , 

and 



||M (n) 1 

|| ?w 

II rp(n) 



(2.55) 
(2.56) 
(2.57) 



HM- 1 !! <^a' d l . 

Combine this with Lemmas 12.21 and 12.31 to obtain the bounds 

||M (nrl || < C 2 5(M (n) )- x = C 2 5(M {0) E tn )- 1 < C 2 C'~ l exp(9t r 
and 

||M< n) || < C^ 1 i(M< n »)-^ = C^ 1 ^MC'E*-)-^ 1 ) 

< Cf^C" - ^ exp[(d - l)0t n ]. (2.58) 
The remaining estimates follow immediately from (|2.49jh ()2.50|) and the equations 

pW=M (n) r 1 "M (0) " 1 , T (n) = M {n) E -5t n ^n-l)- 1 _ ^.59) 

□ 

Proposition 2.5. Choose M = as in (|2.fij) . and suppose a satisfies condition 

()2.29|) . T/ien t/iere zs a constant c-j > stzc/i £/ia£ /or all n G N U {0}, 

"~ -(d-l)V 



1-0 



< |7 (n) | < ci exp[(d - l)9t r 



(2.60) 



C7 exp 
with c\ as in ()2.49|) . 

Proof. The upper bound for |7^"- ) | follows from ()2.49|) . since = (M^)^ and hence 
| 7 W| < ||Af< n >||. 
From ()2.17|) and the Diophantine condition ()2.29j) we have 

| 7 (n) | = 7 exp[(d- l)t n ]\ T p$>a+pM)\ > Cexp[(d-l)t n ]\p^\-^ 1+ ^. (2.61) 



Since 



1 + e 



d-i + e 

(d -!)(!-■ 



(2.62) 
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and 

|j4?l < ll^ll (2-63) 

the proposition follows from the estimate (j2.51j) . □ 

2.8. Hyperbolicity of the transfer matrices. Let 

wf ] = {£ G R d : £ • " (n) = 0} (2.64) 
be the orthogonal complement of the vector 

w (n) = f «^ G R d_ ( 2>65 ) 

Lemma 2.6. For a// £ G u> ( '™~" 1 \ nGN, 

T T (n)-i^ = exp (_ 5 / n ) T (M (n - 1} M ( " rl )£ (2.66) 
Proof. This follows directly from the relation 

E st n T M (n-i)^ = exp( _ 5tn) I* q t j = eX p(-^„) (2.67) 

where £' G ]R d_1 comprises the first d — 1 components of □ 

Proposition 2.7. Choose M = as in f)2.6j) . and suppose a satisfies condition 

(|2.29jl . T/ien t/iere zs a constant A > sucA i/ia£ /or a// £ G a;^™ ^ nGl, 

||T T (.)-^||< l Aexp( _^ n) ||^|| (2 68) 

with 

(p n =(l-0)5t n -d6t n - 1 . (2.69) 

Proof. From Lemma [2.61 

II T T (ri)-1 £|| < exp(-(5t n )||M (n - 1) || ||M^ _1 || ||£||, (2.70) 

and the proposition follows from the bounds ()2.49j) . (|2.5()jl . □ 

Given any positive sequence <p , ipi, . . ., the values t n that solve eq. ()2.69|) with t = 



are 



i=i 

where 

d6 



1 n 

t n = T - § J2( 1 +® n ~ j Vr (2-71) 



/5 = I 3^. (2-72) 



E.g., for constant (p n = if > 0, we have 



t -< U(P {9 = 0) (2 73) 
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2.9. The resonance cone. As we shall see, a crucial step in our renormalization 
scheme is to eliminate all far-from-resonance modes in the Fourier series, i.e., all modes 
labeled by integer vectors outside the cone 



K {n) = {£ G R d : |f -u> (n) | < a n \\£\ 



(2.74) 



for a given a n > 0. 



Lemma 2.8. Choose M = M^ 3 ' as in ()2.fij) . and suppose a. satisfies condition ()2.29|) . 
Then 



II T rW"^i| . 

SU P rprr; < 

-{0} ||f|| 



A 



c 6 o-n-ie' 



q! <5t„ 



e X p[-(l-0)5t n + C ?0t n _ 1 ], (2.75) 



for all fiGN. 

Proof. We write £ = £ x + f 2 ? where 



(n-l) 



(n-l) 



£ 2 e u 



(n-l) 



n) " i mi < ^-lll^ IIIKI 
w (n-l) - n 111 11 IISI 



^(n-l) ||2 

Firstly, 

|| T T^~ :L £ 1 || < || T T^ _1 || ||f 1 || = ||T^ 

since £ G K^-^ and 1 1 <l^_? i> 1 1 = ||(T a (n-i) ) ^ > L Hence in view of 

II T T ( " r ^i|| < c 6 (j n _ 1 exp[( t i - 1)(1 - 6)5t n + d9t n ]\\£\\. 
Secondly, from Proposition 12.71 we infer 

II T T^ -1 f 2 || < ^Aexp [ - (1 - 9)5t n + d9t n _ 1 ] U\\. 

This proves (|2.75j) . 

Remark 2.1. Note that if the t n are chosen as in ()2.73j) . and 



cr n -i < -c 6 1 Aexp(-rf5t n ), 



then 



||T T (n) 1^|| < Aexp(-^)||f| 

for all £ G K {n ~ l \ n G N and <p > 0. 



(2.76) 

(2.77) 

(2.78) 

(2.79) 
□ 

(2.80) 
(2.81) 



3. Renormalization of vector fields 

3.1. Definitions. The transformation of a vector field X on a manifold M by a diffeo- 
morphism if): M — > M is given by the pull-back of X under if): 

i/j*X = {Dif))- l Xoif). 

As the tangent bundle of the d-torus is trivial, TT d ~ T d x M d , we identify the set 
of vector fields on T d with the set of functions from T d to M. d , that can be regarded 
as maps of M d by lifting to the universal cover. We will make use of the analyticity to 
extend to the complex domain, so we will deal with complex analytic functions. We 
will also be considering an extra variable related to a parameter. 
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Remark 3.1. We will be using maps between Banach spaces over C with a notion of 
analyticity stated as follows (cf. e.g. (HI): a map F defined on a domain is analytic 
if it is locally bounded and Gateux differentiate. If it is analytic on a domain, it is 
continuous and Frechet differentiable. Moreover, we have a convergence theorem which 
is going to be used later on. Let {F k } be a sequence of functions analytic and uniformly 
locally bounded on a domain D. If limfc^ +00 F k = F on D, then F is analytic on D. 

Let p, a, b > 0, r = (a, b) and consider the domain 

D p x B r , (3.1) 
where D p = {x G C d : || Imcc|| < p/27r} for the norm \\u\\ = \ui\ on C d , and 

B r = 1 2/ = (yi,...,y d ) e C d : ^2\v*\ < a and |y d | < 6 J . (3.2) 

Take complex analytic functions /: D p x B r — > that are Z d -periodic on the first 
coordinate and on the form of the Fourier series 

f(x,y)=J2h(y)e 2mk ' X - (3-3) 

k&1 d 

Its coefficients are analytic functions fj~: B r — > C d with a continuous extension to the 
closure B r , endowed with the sup-norm: 

\\fk\\ r = SUp ||/fc(l/)||. 

yeB r 

The Banach spaces A p<r and .A' are the subspaces of such functions such that the 
respective norms 

/'■'■ / II J fc|| r 



/ 

p,r 



En/* 

fcez d 

^(l + 27r||fe||)||/ fe || r e*ll 
fcez d 

are finite. Also, write the constant Fourier mode of / G *4p, r through the projection 

E/(l/)= / f(x,y)dx = f (y) (3.4) 



into the projected space denoted by E^4 r . The norm of its derivative D/ is given by 
the operator norm ||-D/ ||r = su P[| s |] r =i \\ D fo9\\r- 

Some of the properties of the above spaces are of easy verification. For instance, given 
any f,ge A' pr we have: 

• \\f{x,y)\\ < \\f\\ p ,r < \\f\\' p , where (x,y) G D p x B r , 

• \\f\\p-S,r < \\f\\p,r With 5 < p. 

In order to setup notations write, according to section|2l =w6l d - {0}, A = 1 
and, for n G N, 



w (n) = 7 (n) X M (n) y J = A n P (n) W = TfcT^W^, (3.5) 

where 

A n = -^e^" 1 ^" and r? n = (3.6) 

7 (n) ln A n _i 1 ^ 
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In the following, we will be interested in equilibria-free vector fields with a "twist" 
along the parameter direction. By rescaling this direction we will find the right param- 
eter which guarantees the conjugacy to a linear flow. For a fixed iigNU {0}, we will 
be studying vector fields of the form 

X(x, y) = X° n (y) + f{x, y), (x, y)eD p x B r , (3.7) 

where / G A p>r and 

^(tf) = « (n) +7 (Brl M (n) »- (3-8) 

(We drop the second coordinate of the vector field because it will always be equal to 
zero - there is no dynamics along the parameter direction.) The linear transformation 
on y deforms the set B r along the directions of the columns of (see (I2.15J) ). In 

particular, its dth column corresponds to u)^ 1 '. 

For the space of the above vector fields we use the same notation A p . r and the same 
norm || • || Pjr without ambiguity. 



3.2. Resonance modes. Given o~ n > we define the far from resonance Fourier modes 
fk as in (|3.3|) with respect to u/ n ) to be the ones whose indices k are in the cone 

J- = {k G Z d : \k ■ uW\ > a n \\k\\}. (3.9) 

Similarly, the resonant modes correspond to the cone 

I+ = Z d -I~. (3.10) 

It is also useful to define the projections 1+ and I~ on A P)T and A' r by restricting 
the Fourier modes to 1+ and I~, respectively. The identity operator is I = 1+ + 1~. 
Moreover, take 

II T^fn+l)- 1 ^!! 

A n = sup — . (3.11) 

fce/+-{0} ll fe ll 

A useful property of the above cones is included in the Lemma below. 
Lemma 3.2. If k 6 I~ and y G B Tn with r n = (a n ,b n ), 

a n <*nQ- b n ^j ^IIIM^II- 1 and b n < i (3.12) 



then 



\k.X° n (y)\>^\\k\\. (3.13) 



Proof. For every y G B Tn and k G I n , 

\ k ■ ( W W + 7 H- 1 M(") ? /)| = |(1 + y d )k ■ W W + 7 (nrl fc • MW(yi, . . . , y d -i, 0)| 

>(l-6„)^||fc||-an|7 (n) l" 1 ||M(")||||fc||. 

Our choice of a n yields fl3.13|) . □ 
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3.3. Basis change, time rescaling and reparametrization. The fundamental step 
of the renormalization is a transformation of the domain of definition of our vector fields. 
This is done by a linear change of basis (coming essentially from the multidimensional 
continued fraction expansion of u> - see section |2J), a linear rescaling of time because the 
orbits take longer to cross the new torus, and a change of variables for the parameter y 
in order to deal with the zero mode of the perturbation. 

Let p n -i, On-ij b n -i > 0, r n _i = (a n _i, 6 n _i) and consider a vector field 

X(x,y)=X n _ 1 (y) + f(x,y), (x,y) G D Pn _ x x B rn _,, (3.15) 

with / G A Pn _ 1)Tn _ 1 . We are interested in the following coordinate and time linear 
changes: 

x^T {nrl x, t^r) n t. (3.16) 

Notice that negative time rescalings are possible, meaning that we are inverting the 
direction of time. In addition to ()3.16|) we will use a transformation on y, a map 
y i— > $> n (X)(y) depending on X in a way to be defined later. 
Therefore, consider the transformation 

L n (x,y) = (^"^^(^(y)), (x,y) G C M , (3.17) 
that determines a vector field in the new coordinates as the image of the map 

X » C n (X) = Vn L* n X. 
That is, for (x,y) G L- l D Pn _ x x B rn _„ 

C n (X)(x,y) = Vn T^{u,^ +y^r 1 M ^<S> n (X)(y) + f o <S> n (X)(y)) 

+ VnT^\f-f )oL n (x,y). 

In order to eliminate the k = mode of the perturbation of X in the new coordinates 

and to normalise the linear term in y to j( n > M^y, using the definitions of and 
rj n we choose 

$ n (X): y^ (ld+ 7 (n - 1) M( n - 1 )~Vo)~^e-^« 2 / 1 ,...,e- rf ^y (i _ 1 ,y rf ), (3.19) 
if possible. Hence, 

C n (X)(x, y) = X° n (x, y) + C n (f - f )(x, y), (3.20) 

where 

C n : f^VnT (n) foL n . (3.21) 
Denote by A M the set of X G -4 p „_ lir „_i such that ||/o||r n -i < I 1 - 
Lemma 3.3. Let r n = (a n , b n ) and [i n -\ > such that 



a n < e° 



dn-l 



bn < &„_! " (l + |7 (n_1) | \\M^ ^ fj, n —l. 



(3.22) 



There exist an analytic map $ n : A Aln _ 1 — > Diff (B r . n , C d ) such that, for each X G A Mn _ 1; 
$ n (X) is given by ()3.19|) and 

$ n (I)(5jc5 r „_ r (3.23) 
In case fo is real-analytic, $ n pf)|K is also real-valued. 
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Proof. For X G A Pn _ u r n -i with ||/ ||r n _i < AV-l and 5 n _i = (5 n -i, 5 n -i) with 

^-i^/Zn-ilT^IHMM -1 !!, (3.24) 
we have by the Cauchy estimate 

\\Df \\ ~ < jl^llrw < 1 

H ■'°llr n _i-« tt -i - § n l " | 7 (n-l)| U^n-l)- 1 !!" 

So, F = Id+7 (n ~ 1) M (n ~ 1) ~ 1 /o is a diffeomorphism on B r _ _j _ Now, if R x < a n _i - 
S n -i ~ fJ-n-u R2 < & n -i ~ #n-i - A*n-i and R = (R 1 ,R 2 ), we have B R C F^^^J 
and F~ 1 {Br) C B r j . Therefore, $ n (X) as given by ()3.19|) is a diffeomorphism on 

S rn by choosing i? = (e~ d5tn a n , b n ), and thus we get ()3.23|) . In addition, X i— >■ $ n (X) is 
analytic from its dependence on / . When restricted to a real domain for a real- analytic 
/o, $ n (X) is also real-analytic. □ 

Let the translation R z on C 2d be defined for z G C d and given by 

i-> (05 + 2,2/). (3.25) 

Notice that we have the following "commutative" relation: 

L* n R* z = R* T(n)z L* n , z G C d . (3.26) 

This also follows from the fact that $ n is unchanged by the introduction of the trans- 
lation R z . 

3.4. Analyticity improvement. 
Lemma 3.4. If 5 > and 



P ' n <f^-5, (3.27) 



1 n— 1 



then C n as a map from (IL^_ 1 — E)^4 /9n _ 1)7 . n _ 1 fl A Atn _ 1 znto (I — ~E)A' p , rn continuous 
and compact with 



\C n \\<\vn\\\T^\\ (l + y)- 



(3.28) 



Remark 3.5. This result means that every vector field in I^_ 1 ^4. Pn _ li r n _i ^ ^^-n i- e - 
a function on .D Pn _ 1 x S r . n _ 1 into C d , has an analytic extension to 1 D p ' x B Tn l . 

Proof. Let / G (I+- X - E)^,^ n A^. Then, 
Wf° L n\\'^r n < £ (l + 27r|| T TW- l fc||) ||/ fc oM n || r „e^~^)ll T ^ rlfc ll. (3.29) 

By using the relation £e _<5 £ < with £ > 0, ()3.11|) and (|3.23|) . we get 

ll/oi»li; ; , r „<(l + 2^) £ INk-x^" 1 ^ 11 * 11 

d-M (3.30) 

<(1 + 2tt/5) ll/IU.^. 

Finally, < M ||rM|| ||/ o Lj,, ^ 

The above for D p i n x _B rn is also valid for x B rn , Q > p' n but satisfying a similar 
inequality to ()3.27|) . Therefore, C n = To J , where J ': — E)^4 pn _ ljrn _ 1 — > *4.£ j7 . n is 

bounded as C n , and the inclusion map X: >A£ jrfi — > .A p / rn is compact. □ 
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For < p" n < p' n , consider the inclusion 

1 n '. A n l „ ► A n ll r (3.31) 

by restricting X G A', _ to the smaller domain D a n x B Tm . When restricted to non- 
constant modes, its norm can be estimated as follows. 

Lemma 3.6. If<p n > 1 and 

O<p:< P ;-log(0 n ), (3.32) 

then 

||X„(I-E)|| <<P~\ (3.33) 
Proof. For / G (I — E)A' pl >rn , we have 

fc^O 

□ 

3.5. Elimination of far from resonance modes. The theorem below (to be proven 
in Section fA. 1)) states the existence of a nonlinear change of coordinates U, isotopic to 
the identity, that cancels the I~ modes of any X as in ()3.7|) with sufficiently small /. We 
are eliminating only the far from resonance modes, this way avoiding the complications 
usually related to small divisors. We remark that the "parameter" direction y is not 
affected by this change of coordinates. 

For given p n ,r n ,e, v > 0, denote by V £ the open ball in A pn+V ^ n centred at X° with 
radius e. 

Theorem 3.7. Let r n be as in dSH, o n < \\u>^\\ and 

° n min \ TZ> ^TTTr^n t ■ ( 3 - 35 ) 



(3.36) 



42 [4tt' 72||g>) 
For all X G V £n there exists an isotopy 

U t : D Pn x B Tn -> D Pn+u x B rn , 

(x,y) !-»• (x + u t (x,y),y 
of analytic diffeomorphisms with u t in A Pn Tn , t G [0, 1], satisfying 

I-U t *X = (l-t)I~X, U = ld. (3.37) 

This defines the maps 



and 



X ^ Id +w t 



W t : V £n - © (1 " 

X ^ u*x 



which are analytic, and satisfy the inequalities 

42f 

lin t (x) - id n;„ >rn <— \\i-f\\ Pn , rn 

On 

\\U t (X) - X°j pn ,r n <(3-t)\\f\\> pM . 

IfX is real- analytic, then < d t (X){R 2d ) C M M . 



(3.38) 



(3.39) 



(3.40) 
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Remark 3.8. Further on we will be using the above result for t = 1. So that all far 
from resonance modes are eliminated. 



Recall the definition of the translation R z in Ij3.25|) . 

Lemma 3.9. In the conditions of Theorem \3.7\ if x &R d and X G V £n , then 

il t (X o R x ) = R~ l o ii t (X) o R x (3.41) 

on V Pn>Tn . 

Proof. Notice that R x (D Pn x B Tn ) = D Pn x B Tn . If U t = iit(X) is a solution of the 
homotopy equation (|3.37j) on D Pn x B rn , then U t = R~ l oil t (X) o R x solves the same 
equation for X = X o R Xy i.e. I~X o U t = (1 — t)I~X, on D Pn x B Tn . □ 

3.6. Trivial limit of renormalization. Let a sequence of "widths" < o~ n < 1 of the 

resonance cones 1+ be given. The nth step renormalization operator is thus 

Tl n = U n o J n o C n o 7^„_i and 7£ = Ho, 

where U n is the full elimination of the modes in I~ as in Theorem 13 . 71 (for t — 1). Notice 
that 7£ n (X° + t>) = for every v G C d . From the previous sections the map Tl n on 
its domain is analytic. Also, in case a vector field X is real-analytic, the same is true 
for Tl n {X). 

Fix the constants v and 5 as in Theorem 13 . 71 and Lemma HOJ respectively, and choose 
< A < 1. Take 

9 n = min e n , " 1 , A" ^ g^g^ 

by assuming that the sequence of times t n guarantees that n > 0. Now, write 

n 

B n = Y[A (3.43) 

t=0 

with Ai given by ()3.11|) . By recalling the inequalities ()3.27|) and ()3.32|) we choose, for a 
given p > 0, 



5 

n— 1 n— 1 



where 



i=0 i=0 



(3.44) 



max |2|r/ n | ||T( n )||(l + 27rcT 1 )^pi, l| > 1 (3.45) 



is to be used in Lemma ETol 

Define the following function for every u> G M. d associated to the choice of o~ n : 

+oo +oo 

= X A log + (5 + i/) B i- ( 3 - 46 ) 

i=0 i=0 

The convergence of the renormalization scheme now follows directly from our con- 
struction. 

Theorem 3.10. Suppose that 

B(u) < +oo (3.47) 
and p > B(u>) + v . There is K,b > and r n = (a n , b n ) with a n > and b n > b > 0, 
such that if X is in a sufficiently small open ball around X° in A p , ro , then 
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(i) X is in the domain oflZ n and 

\\n n (X) - n n (X°)\\ Pntrn < KQ n \\X - X°|| p , ro , n G N U {0}, (3.48) 

(ii) for each \s\ < b there exists in B r . n l C C d the limits 

p s n (X)= lim $ n (^„_ 1 (X))...$ m (^ m _ 1 (X))(0,...,0, S ) (3.49) 

m— »+oo 

and 

lim |b*(X)-(0,...,0, S )|| =0, (3.50) 

n— >+oo 

(iii) the map X i— > p n (X) is analytic and takes any real-analytic X into an analytic 
curve s t— > p s n {X) in M. d . 

Proof. Let £ > and po = p — v — £ > such that po > ^(<*0- Hence, by ()3.44|) . we 
have R > satisfying p n > RB~\ for all n G N. 

Denote by c the radius of an open ball in A Ptro centred at X° and containing X. If 
c < £o we can use Theorem 13.71 to obtain TZo(X) G Iq *4 POiro with r = (ao, &o) satisfying 
flSH) and 

l^o(X) - ^ (X°)|| po , r , < 2||X - X°|| p+?iro < 2r X ||X - X°|| p , ro . 

Let if = 2(£6 )~ 1 and assume that c < K _1 min{6 (l - A), ~ - 6 }- So, (j3~3%jl holds 
for n = 0. 

Now, with n G N we choose the following r n : 

A \ M n )| ^ . 

On = On ^ ~ &0 J ||Af( w )|| aIld & " = & ° _ CK ^' ^' 51 ' ) 

so that 1/2 > b n > b = bo — cK(l — A) -1 . The inequalities in (j3.12|) follow immediately. 
Moreover, ()3.22|) is also satisfied with p n -\ = cKQ n _i because 

d5tn " >(-- 6 )(1 + |7 (n_1) | HM^- 1 ^ 1 !^,,.! 



a n _ x - e "a H _ , o 



> (1 + 17^-^1 \\M {n - 1] 1 \\)cKQ 



n—li 



b„ , - b n = cK\ n - 1 



(3.52) 



Suppose that X n _i = TZ n -i(X) G I n "_i^lp n „ 1 ,r„_i and 

||X n _i - -X'n-ill/^-i.j-n-i ^ X6„_i||X - X°|| Pir . 

Since ()3.22|) holds, Lemmas 13 . 31 and I3~4l are valid and, together with ()3.20|) and Lemma 
13. 6[ can be used to estimate X n o £„(X n _i): 

|| I n o C n (X n ^) - X%„ rn < \ Vn \ \\TW\\(1 + 2n5- l )r n l KQ ri ^\\X - X°|| p , ro 

1 (3-53) 

= -X6 n ||X — X || Pjro . 

This vector field is inside the domain of U n as (j3.12|) and \cKQ n < e n are satisfied. 
Thus (EHHj) follows from (COnjl . 

Denote by /q the constant mode of the perturbation term of X n . By Lemma 13.31 
$ n (X„_x) : B rn -> 5 rn _ 1 is given by 

y h-> (Id +<7 n ) diag(e~^, . . . , e"^, 
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where 

9n = (id+^M^-'ft^y 1 - Id (3.54) 
is defined on B r r with r' n = (e~ d5t "a n , b n ). So, for z G B r > n there is £ G B r / n such that 

g n (z) = [/ + 7 ( "- 1) M("- 1 )" 1 D/ ( "- 1) (0]- 1 [^-7 (n - 1) M("- 1 )"Vo (n_1) (0)] -z 



.[/ + 7 (n~i) M (n-i) 1 D/i"" 1) (0] 4 7 ( " 4) M M * + /o (ri_1) (0)] 

(3.55) 



and 



< ^' 1)|l|M( 7 irl|l ( a (\y n \\ W D ft%> + ll/o^ilA 



(3.56) 



The choice of r n means that 



minK-x - e~ d5t "a n , b n ^ - b n } > min [ ^(l-")^ ^ } ■ (3-57) 

By using ()3.42|) and the Cauchy estimate, 

II f( n_1 )|| \n-l 

IP/o n_1) llr' < — T " J ° J, rr « i- (3.58) 

jo "^-minK^-e-^o^ti-y | T (n-i)| j| tv^C— i) " x || 

Thus, 

||^n|k<A n - 1 . (3.59) 
Writing y s = (0, . . . , 0, s), by induction we have 

m 

$ n (X„_i) . . . $ m (X m „ 1 )(?/J = V. + Yl diag(e-^- 1 -*- 1 ), . . . , e^-"'—), l)^(^), 

i=n 

(3.60) 

for some £ k G B r < . Therefore, from ()3.59|) . there exists p s n {X) G C d unless X is real 
which clearly gives p s n (X) G IR d . In addition, 

\\ P s n {x) - y s \\ < WmWri « y=x- ( 3 - 61 ) 



The maps X \— > are analytic since the convergence is uniform. Lemma 

gives us the nested sequence <& n (X n _i)(B rn ) C B r . n l . So, as y s G n ie N-B ri , it follows 

Remark 3.11. The above can be generalised for a small analyticity radius p by consider- 
ing a sufficiently large N and applying the above theorem to X = IA^Cn ■ ■ ■ U\C-JAq{X), 
where X is close enough to X°. We recover the large strip case since is of the order 
of -B/v^!- It remains to check that p^ > B(lj^) + v. This follows from the fact that 
B(«W) = B N \\B{u) - B N (u)] where B n (uj) is the sum of the first N terms of 
so that i3jv(uj) — > iB(tt>) as TV — ► +oo. 

Lemma 3.12. If uj = ( ") m R d is diophantine, i.e. a satisfies (J2.29)) wift exponent e 
(Wated to 6> fry P^Hjl and to /3 fry (pT72jl j ; i/ien (j3~4Tjl verified. 
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Proof. Let us set 8t n = £t„_i, o~ n = exp(— c<5t n ), n > 1, where positive constants £, c 
will be chosen later. Obviously, t n = (l + £)£ n _i = [(1 + an d St n = (l + £)<5t n _i. 

We shall assume that 

c< + (3.62) 
so that cr„_i exp(dd~t n ) = exp(— c5t n _i + d5t n ) = exp[(d — c/(l + £))6t n ) is much larger 
than A given by Proposition 12.71 Hence, using ()2.75|) we have 

We next estimate ||^ (n) || and e n . It follows from (JHIEJ) that ||u/ n )|| < ||M( n )|| |7 (n)_1 |. 
Thus, using (l2~o7IJ) we have 



< exp 



5t r 



(3.63) 



|o; (n) || < exp 



exp 



l + £ 



1-0 £ 



5t T 



(3.64) 



Since ||w ( ' n ' ) || > 1 one gets from ()3.35|) that e n ~ cr^/||a>^|| which together with ()3.64|1 
implies 



exp 



-2c 



8 



5t r 



<C e n <§; exp(— 2c5t n ). 



(3.65) 



1-0 e 

Here X ~ Y means that there exist two positive constants C\, C% > such that C{Y < 
X < C 2 Y. Using again (J23HJ), (M) we get 



Cn|7 



(n)| 



|| 

Also, since ||M( n+1 )|| > ^ n+1 \ 



3> exp 



-c — 



fll + C 



1-9 £ 



Str 



< exp[-(c + d)(l + f)<**r 



We shall assume that c and £ are chosen in such a way that 



1-0 e 



>-(c + d)(l + 0, 



(3.66) 

(3.67) 

(3.68) 
(3.69) 

i-e e- { " (370) 

Finally, we want A n to be small and, hence, require the exponent in ()3.63|) to be 
negative 

(3.71) 



so that 



(n)\ 



(n+l)| 



> 



0"n|7 



Wl 



||M( n )|| ||M( n+1 )||exp(dt n+ i) ||M(™)||' 
Inequality (|3.68jl is equivalent to the following condition 



c > 



e 1 + ^. i+{ 



d. 



c d9 

— + d-(i-e) + j<o. 



Suppose that conditions (|3.62jl . (|3.7()jl . ()3.71|) are satisfied. It follows immediately 
from the estimates above and (jSH), flZHS) , HI , (fi^fy . (|2~Knj) that 

I log e n | , | log e n _! | , log ||r(") || , | log |^ n | | < st n . (3.72) 



At the same time 



B n = J[Ai < C"exp(-at n+1 ) 



(3.73) 



i=0 
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where 

a = -^--d+(l-9)-^->0. (3.74) 

Since _B n decays exponentially with t n and log <p n grows at most linearly the series ()3.47|) 
converges. To finish the proof it is enough to show that conditions ()3.62j) . ()3.7()|) . ()3.71j) 
can be satisfied. Indeed, since < 9 < 1 we can choose £ so large that 1 — 6> — d9/£ > 
and 

x _g ^ 2 <* --7— rf<0 - ( 3 - 75 ) 

It is easy to see that all three inequalities (|3.62j) . ()3.70|) . ()3.71|) are satisfied if + — 
/?) < c < (1 + f )d, where = 1-6- d6/£ > 0. □ 

3.7. Analytic conjugacy to linear flow. As a consequence of Theorem I3.1(J[ we 
obtain an analytic conjugacy between a vector field and the linear flow, thus proving 
Theorem ll.il In the following we always assume to be in the conditions of Section I3~6l 
Let r = r and 

A = {X G A p y. ||A-A°|| Ar < c} (3.76) 

inside the domain of lZ n for all n G N U {0}. By taking A G A, we denote X n = 
K n {X) G I+Aw„ so that 

X n = \ n (U o L x o U t ■ ■ ■ L n o U n )*(X), (3.77) 

where C4 = llfc(X fc £ fc (A fc _ 1 )) is given by Theorem 13 .71 for t = 1 at the fcth step and L fc 
is the linear rescaling as in ()3.17|) for A*._i. 
Denote by the coordinate change 

V n . (x,y) i- (p('^ 1 a; ,$ 1 (A )...<l> n (A„_ 1 )( 2 /)) (3.78) 

and set Vo = Id. Thus, L n = V~_ x o V n and 

X n = A n (K o f/ n )*(K-i o f/ n „! o • • • (V x o U x o V^yU^X). (3.79) 

In particular, the y-coordinate is only transformed by the second component of V n . 
Notice that if X n = A° for some iigN, 

y = $!(A ) . . . $ n (X n _ 1 )(0, . . . , 0, s) G C d , 

with \s\ < b, corresponds to the parameter for which X is conjugated to (1 + s) u}^ 71 '. 
The parameter value for the general case X n — A° — > as n — > +oo is p s {X) = p\[X). 

Lemma 3.13. There is an open ball B about A in A such that we can find a sequence 
R n > satisfying R_i = p, 

R n + 2vr42Aey 2 ||A - X°\\ PjT < R n ^ < j^^, X G B, (3.80) 

and 

lim R-'Ql/ 2 = 0. (3.81) 

n— >+oo 

Proof. Let = minp n . It is enough to check that Ql/ 2 <C A n p* niLi ll^^ll -1 with 

1/2 

< A < 1 and taking R n = c\~ n Q n for some positive constant c. This immediately 
implies (13.81)) and ()3.8U|) by considering a small enough upper bound for \\H—H \\ p ^. □ 
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Let Diff per ( J D c ,C d ), ( > 0, be the Banach space of Z d -periodic diffeomorphisms 
g: Dq — > C d with finite norm \\g\\(; = J^ fc ||5'fc|| e ''" fc ") where g^ G C d are the coefficients 
of the Fourier representation. It is simple to check that \\g o P^\\ Rn < \\g\\ P „- 

Denote by u n the analytic function 

u n . A ^Diff per (D Pn , C d ) 

Xh->il n (2 ft A l (X n _ 1 ))(.,p- +1 (X)). 

As p s n+1 {X) G B rn , D Pn x {y = p s n+1 (X)} is inside the domain V Pn)Tn of il n (X n £ n (X n _i)) 
given in Theorem 13.71 Now, for each X, define the isotopic to the identity diffeomor- 
phism 

W n {X) = pW 1 o Un {X) o ?W, (3.83) 

on P^ 1 D Pn . If X is real- analytic, then W n (X){R d ) C IR d , since this property holds 
for Un{X). We also have W n (X°) = Id. 

Lemma 3.14. For a// n G N U {0}, W n : B — ► Diff per (D Rn , C d ) zs analytic satisfying 
WjXy.D^^D^ and 

\\W n {X) - Id Ik < 42K6y 2 ||X - X°|| p , n X G £?. (3.84) 

Proof. For any X G A, in view of ()3.40|) we get 

\\w n (x) - id \\ Rn = o [ Mn (x) - id] o pwy^ 

We can bound the above by ()3.84|) . 
Now, for x G D Rn and X G S C A, 

|| \mW n {X)(x)\\ < || lm(W n (X)(x) - x)\\ + || lmx\\ 

< \\W n (X) - Id Ik + i4/27T < p n _ 1 /2 7 r. 

So we have W„(X) : D Rn — > -D/ ?n _ 1 and W„(X) G Diff per .(Di ?n , C d ). From the properties 
of iln, W n is analytic as a map from 5 into Di& per (D Rn , C d ). □ 

Consider the analytic map H n : B — > Diff per (D/j n , C d ) defined by the coordinate trans- 
formation H n (X) : Z}# n — > D po as 

#„(X) = W (X) o ■ ■ ■ o W n (X). (3.85) 

Lemma 3.15. There exists c > snc/i t/iat /or X G -B and n6N, 

||i/„(X)-if„_ 1 (X)||^ <cey 2 ||X-X°|| Pir . (3.86) 

Proof. For each fc = 0,...,n — 1, consider the transformations 

G k (z,X) ={W k (X)-Id) o (Id+G fe+1 (z,X)) + G k+1 (z,X), 
G n (z,X) =z(W n (X)-ld), 

with (z, X) G {-2 G C: \z\ < 1 + d n } x 5, where we have d > such that 

c' 

ey 2 ||x-xo|| Pir 

If the image of D Rn under Id +Gk+i(z, X) is inside the domain of Wk(X), or simply 

\\G k+1 (z,X)\\ Rn < (R k -R n )/2n, 
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then G k is well-defined as an analytic map into Diff per (Z}^ n , C d ), and 

\\G k (z,X)\\ Rn < \\W k (X)-ld\\ Rk + \\G k+1 (z,X)\\ Rn . 

An inductive scheme shows that 

\\G n {z,X)\\ Rn <(R n ^-R n )/2n, 
n-l 

\\G k (z,X)\\ Rn < \\Wi{X) - Id |U + |z| \\W n {X) - Id |k 

i=fc 

<(R k -X - Rn)/27T. 

By Cauchy's formula 

\\H n (X) - H n ^{X)\\ Rn = \\G (l,X)-G Q (0,X)\\ Rn 

1 r G (z,X) 



2m 



dz 



z\=l+d n /2 Z \ Z 



1 

R n 



and 



\H n (X) -if^pOlk < A sup ||C7 (^,A)|| Rn 

U n | 2 |=l+d n /2 



«ey 2 ||x-x°[, 

□ 



Consider <^ er (M d , C d ) to be the Banach space of the Z d -periodic C 1 functions between 
i d and C d with norm 

C i =maxmax||D*/(aj)||. (3.87) 

fc<i 



Lemma 3.16. There exists C > 0, an open ball B' C B about X° and an analytic map 
H: B' -> DifT per (R d , C d ) suc/i i/mt /or X G B' , H(X) = lim n ^ +00 # n (A) and 

\\H(X) - Id || c i < C\\X - X°\\ p>r . (3.88) 

IfXeB' is real- analytic, then H(X) G Diff per .(M d , R d ). 

Proof. As the domains D Rn are shrinking, we consider the restrictions of W n (X) and 
H n (X) to M d , and estimate their C 1 norms from the respective norms in Di& per (D Rn , C d ). 
More precisely, for any X G B, making use of Lemma 13. 15| 

\\H n (X) - H n ^(X)\y <max sup \\D k [H n (X)(x) - H n ^(X)(x)}\\ 



k<l 



xeD 



Rn/2 



(3.89) 

<— \\H n (X) - H n ^(X)\\ Rn , 

which goes to zero by (j3.81)l . Notice that here we have used Cauchy's estimate H-D 1 ^^ < 
(27r/e5)||0|| c+ 5 with C, 5 > 0. 

Therefore, it is shown the existence of the limit H n (X) — > H(X) as n — * +oo, in the 
Banach space Cp er (R d , C d ). Moreover, ||if(A) — Id H^i <C \\X — X°\\ PiJ .. The convergence 
of H n is uniform in B so H is analytic. As the space of close to identity diffeomorphisms 
is closed for the C 1 norm, H(X) is a diffeomorphism for any X sufficiently close to A , 
i.e. A G B' . The fact that, for real-analytic A, H(X) takes real values for real 
arguments, follows from the same property of each W n (X). □ 

To simplify notation, write TT y X = X(-,y). 
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Lemma 3.17. For every real-analytic X G B' and \s\ < b, [H (X)]* (tt p s ( X )X) = (1+s) U3 



on iii d . 



Proof. For each n£N the definition of H n (X) and ()3.79j) imply that 

H n (X)*{-K pS{x) X) = A; 1 7r pS{x) V;- 1 *(X n ). (3.90) 
The r.h.s. can be written as 

A; 1 pW" 1 [a;W+ 7 W- 1 MW$ n (X n _ 1 )- 1 • ■ ■ $i(X )" V (X)]+A; VpoKT^-X") = 
= (1 + s) u + X-'P^^M^p^X) -suj + A„- 1 V(X) V;- 1 *(X„ - X°). (3.91) 
Its terms can be estimated, for x G R d , by 

llA^vwKT 1 *^-^)^)!! < Ia; 1 ! wp^w \\x n -x° n \\ Pn ,r n « e^ 2 , (3.92) 

and 

_^|| P H- M (n)^ +i(X) _ (0, . . . , 0, S )]|| = J^^\m P S n+1 (X) ~ (0, . . . , 0, S)] || 

«||^ +1 (X)-(0,...,0, S )|| 

(3.93) 

which is controlled by (|3.5()j) . Consequently, the limit of ()3.9()|) as n — > +00 is (1 + s) U3. 
Using the convergence of H n we complete the proof. □ 

Lemma 3.18. If X e B' and x G R d , then 

H(X o R x ) = R' 1 o H(X) o R x , (3.94) 

where R x : z 1— > z + x is a translation on C d . 

Proof. For each n G N, (I3~4T1) and (jQBj) yield that W n (X o i^) = U n {X) o and 
C n (X o i?^) = £ n (X) o R T(n)x . This implies immediately that 

n n (X o i^) = TZ n (X) o i2 p(n)B . (3.95) 

Next, from a simple adaptation of (j3.41|) and the formula R P (n) z = P^R Z P^ 1 for 
z G C d , we get 

W n (X o R x ) =P^ ott an (C n Kn-i(X O i^)) O ?W 

^ ^ (3.9oJ 

= J R; 1 oiy„(X)oi? £C . 

Thus, H n (X o i?^) = i?" 1 o H n (X) o i? x . The convergence of If n implies (13.94)) . □ 

Theorem 3.19. Ifv G Vect^(T d ) is sufficiently close to uj, then there exists an analytic 
curve s 1— > p s G M. d for \s\ < b, and h G Diff a; (T d ) homotopic to the identity such that 

h*( v +p s ) = (l + s)u>. (3.97) 

The maps v 1— ► p and v *—>■ h are analytic. 

Proof. The lift v to R d of v is assumed to have an analytic extension in D p . Consider 
the real-analytic vector field 

Y(x,y) = v(x) + y 

in Ap, r - Suppose that v is close enough to u> such that Y G B' and Y oR z g B' for some 
77 > and 2 G .D^. Then, the parameter p s = p s (Y) G R d and the C 1 -diffeomorphism 
h = H(Y) mod 1 verify ((337)1 . 
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We now want to extend h analytically to a complex neighbourhood of its domain. 
Take h(z) = z + H(Y o R z )(0), z G D v . The maps z ^ Y o R z and X i-> F(Jf ) are 
analytic and Cp er (M d , C d ) 3 g \— > g(0) is bounded. As /i involves their composition, it is 
analytic on the domain D v and Z d -periodic. From ([3.940 . for any x G R d , we have 

h(x) mod 1 = (a; + ^ o #(y) o ^(0)) mod 1 

= (x + H(Y)(x) - x) modi (3.98) 
= h(x). 

The extension of h is a complex analytic diffeomorphism, thus h is a real-analytic 
diffeomorphism. □ 



4. RENORMALIZATION OF HAMILTONIAN FLOWS 

4.1. Preliminaries. Consider the symplectic manifold T*T d with respect to the canon- 
ical symplectic form X^=i dy% A dx{. As the cotangent bundle of is trivial, T*T d ~ 
T d x R d , we identify functions on T*T d with functions on T d x R d . By lifting to the 
universal cover, we consider functions from ~R 2d into R and extend them to the complex 
domain. 

Let Q be a neighbourhood of R d x {0} in C 2d . A Hamiltonian is a complex analytic 
function H: Q — ► C, Z d -periodic on the first coordinate, written on the form of a 
Taylor- Fourier series 

H(x,y) = J2 H k ,„y»e 2 * ik x , (x, y) G O, (4.1) 
(fe,i/)e/ 

where I = Z d x (NU {0}) d , if fe)1/ G C and y" = y^ 1 . . . 

Let the positive real numbers p and r be given in order to determine the domain 



where 



V Ptr = D p x B r , (4.2) 

D a = \x G C d : II Imajll < p/2tt} and 
pi II (4.3) 

Br = {yeC d : ||l/||<r}, 

for the norm ||u|| = J^ =1 on ( ^- a! - Consider the Banach space A p>r of Hamiltonians 
defined on Q = T> PiT , which extend continuously to the boundary and with finite norm 

\\H\\ p , r = J2 \HuAr M] ^ M - (4.4) 
(k,u)ei 

Similarly, take a norm on the product space A 2d r = A p>r x ■ • • x A p>r given by || (Hi , . . . , H 2 d) \\ P)r 

Ym=i \\Hi\\p,r- Using this we have the Banach space A' of the Hamiltonians H G A PlT 
with finite norm 

ll-f llp,r = ll-f ll/a,r + II VH\\ Pir . 

A useful property is the Cauchy estimate: for any 5 > we have 

1 1 Si-Hi Ur < ~rl|-H||p+(5,r, H G A p+ s, r , 1 < i < d, 

/ ' (4-5) 

H^^ILr < ^||-ff||p, r+ tf, H E A P: r+5, d+l<j<2d, 
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where dk denotes the partial derivative with respect to the kth argument. In particular 

\\H\\' P> r < U + W H \\p + S,r + 6. (4.6) 

The constant Fourier modes are written by the projections 
E F= [ F(x,0)dx = F 0fi and EF(y) = [ F(x,y)dx = Vf / '• (4-7) 

The space where KF lies is denoted by KA r and the natural induced norm is || ■ || r . 
Similarly, we define KA' r with norm || ■ \\' r . Recall Remark 13.11 about the notion of 
analyticity of maps over the above spaces. 

In the next sections we will be studying Hamiltonians close enough to H° as given 
by ()1.7j) for a given choice of uj e M. d — {0} and an invertible symmetric d x d matrix 
Q. The norm of matrices is given by ||Q|| = maXj = i.. >( j X)i=i \Qi,j\-> w h ere Qi,j are the 
entries of the d x d matrix Q. 

4.2. Change of basis and rescaling. The following transformations leave invariant 
the dynamics of the flow generated by the Hamiltonian, producing an equivalent system. 
They consist of 

• an affine linear symplectic transformation of the phase space, 

L n . (x,y) i ► (TW _1 aj, T T^y + b n ), {x,y) e C M , (4.8) 

for some b n e C d , 

• a linear time (energy) change, 

H i ► 77 n i?, (4.9) 

where r\ n is as in (|3.6j) . 

• a linear action rescaling, 

H^—H(-,ii n -) (4.10) 

with a choice of \i n > to be specified later on, 

• and the (trivial) elimination of the constant term 

H^(I-E )H. (4.11) 

Recall the definition of the sequence lj^ given by (|3.5jl . For neN and p n _i,r > 0, 
consider Hamiltonians of the form 

H(x, y) = H Q n _ x {y) + F(x, y), (x, y) e V Pn _ 1>r , (4.12) 

where 

Hl_ x {y) = ■y+ 1 - T yQ n - iy (4.13) 

and H$ = H°, with Q n -\ being a d x d symmetric matrix. By performing the transfor- 
mations above, we get a new Hamiltonian that is the image of the map 

H i ► C n (H) = (I — E )^iJ o L n (-, /i n -). 

In order to simplify notation, we write 

§ n {y) = Li n T T^y + b n . (4.14) 



(I-Eo)^ 

Mr, 



■ $ n (y) + l - T $ n (y)Q^<!> n {y) + F o L n (x, fi n y) 
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So, for any (aj, y) E L^ x V Pn _^ r , 
C n (H)(x,y) = 

f ' ' 1 ■ (4-15) 

By decomposing F — (I — K)F + F and using the Taylor expansion of F : 

F o$ n (y) =F (b n )+^ n T VF (M T TWy + ^ T yTWD 2 F (bO T r (n) y + T n (^), (4.16) 

withT„(y)=0(||y|| 3 ),weget 

£>n{H){x, y) =u>W -y + Vn [ T 6 n <3„-i + T VF (b n )] T T^y 

+ ^ T yT<"» [Q„_i + D 2 F (6 n )] T T^y (41?) 

+ ^T n (y) + ^(1 - E)F o // n y). 

To "normalize" the linear terms in y of E£ n (i7) equal to ■ y, we choose b n inside 
the domain of VF such that 

Qn^K + VF (6„) = 0. (4.18) 

Hence, 

C n (H)(x,y) =H°(x,y) + £ n (F )(y) + £ n (F-F )(x,y). (4.19) 
Here we have introduced the operators 



and 



C n :F ^^T n (4.20) 

fin 

l n . (I-E)F^^(I-E)FoL n (. )/v ) (4.21) 

defined in KA r and (I — E,)A Pn _ 1 , r , respectively. Moreover, we have used the symmetric 
d x d matrix 

Qn = Vn^nT {n) [Q„_l + D 2 F (b n )] T T^. (4.22) 
Denote by A 7 the set of H e A Pn _ 1>r such that ||-Fo|| Pn _ li7 . < 7- 
Lemma 4.1. If det(Q n _i) 7^ and 

2 

7n = T , , (4.23) 
t/iere is b n G C 1 (A 7n ,C d ) suc/i that, for all H e A 7n , 6 n = b n (H) satisfies ()4.18j) and 

||M^)ll<(2A)ll^-illPo||,<^ (4.24) 
Moreover, det(Q„) 7^ and 

, ||T (n)_1 || || T T (n)_1 || 

"^"-MnkKltelll^-Sll^ollr)" (425) 

In t/ie case Fq real- analytic and Q n -i is real, then b n (H) e M d and Q n zs a/so reaZ. 
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Proof. Consider the differentiable function !F(H, b) = b + Q~_i7F (b) defined on 
A 7n x B r / 2 . Notice that jF(iJ°_ 1; 0) = and the derivative with respect to the sec- 
ond argument, 



D 2 F(H, b)=I + Q-\D 2 F {b), (H, b) e A ln x B r/2 , 



admits a bounded inverse because 

||-D 2 -Fo||r/2 = max ||9,-VF || r /2 

d+l<j<2d J ' 

<(4/r)||VF || 3r/ 4 (4 26) 

< (16/r 2 )||F || r 

< iiQ^iir 1 

by the Cauchy estimate. The implicit function theorem implies the existence of a C 1 
function b: H i— > b(H) in a neighbourhood of H^_ l such that 

T{H, b(H)) = b(H) + Q-^VF (b(H)) = 0, 

i.e. a solution of 1)4.18]) . Notice that for any H G A 7n the operator Id— J-"(H, •) is a 
contraction with a unique fixed point b. Hence the domain of the C l function H \— > 
b(H) = b is extendable to A 7n and thus ()4.24|) . Assuming F to be real-analytic and 
Q n -i with real entries, the same argument is still valid when considering B r / 2 D M. d . 
Thus, b(H) is real and Q n is a real symmetric matrix. 
From (jOHjl . 

||q;V^o(M#))II < i, ^eA 7n . 

Hence, A = <3 n _i[I + Q^i-D 2 i 7 o(&n(-ff))] is invertible. Moreover, 

II^H^l/dlQ^II^-II^Foll^). (4.27) 

Now, Q- 1 = (rinfin)- 1 T T^~ l A- l T^~ X , thus (j4~25j) . □ 



Lemma 4.2. If r < r' and 

then C n : E^4 r fl A 7n — > E^4y and 



Mn < — ,, -r m , mi ) (4.28) 



|£ II < u 2 |r7 | f 1 + 1-] ( 4/ ll TT(n) ID 3 (4 29) 



Proo/. Let # G A 7n , i? = 4r , Mn|| r TT(n)|| > 1 and the map 

/: {z G C: \z\ <R}^C d 

z^F (zfi n T T^y + b n (H)). 
Hence T n as in (|4.1fij) can be written as 

/(l) - /(0) - DM - itf/(0) = ^ ^ jjg^fa. 



(4.30) 



Therefore, 



I "V^ II' 

Z7T 



1 '• 



\ Z \=R Z 3 (Z - 1) 



dz 



^ mtp u sup ll*o(^ T T (n) • +bn(H))\\' rl . 
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For \\y\\ < r', in view of f)4.24j) . 

sup \\zfi n T T^y + b n (H)\\ <RHn\\ T T {n) \\r'+\\b n (H)\\ < r/2, 

\z\=R 



and 



sup \\F ( Zf x n T T^ ■ +b n (H))\\' r , < \\F \\ r/2 + Rp, n || T T^|| ||VF || r/2 

\z\=R 

n n 1 .. ., / 1 \ „ „ ( 431 > 

< II -Poll r/2 + 7^7ll^0||r < I 1 + — I ll^b||r- 



Thus, ||T n ||;, < (1 + l/2r')[R 2 (R - l^H-Follr and 

im^o)|i;,^iitx<M( 1 + _L) 



(4rVn| 


T T (n) | ) 3 


^•2 


4r' 


|/^n 


|| W||) 



Or- 



□ 

Finally, we have the relation for each z G C d : 

R z ° L n = L n o R T (n) z , (4.32) 

because EH o R z = EH. 

4.3. Far from resonance modes. Given o~ n , r n > 0, we call far from resonance modes 
with respect to u/ n ) the Taylor- Fourier modes with indices in 

J- = {(fe,i/) G I: |w (n) ■ fc| > er n ||Jfe||, \\u\\ < r n ||fe||}. (4.33) 

The resonant modes are in = I — I~. We also have the projections 1+ and I~ 
over the spaces of Hamiltonians by restricting the Taylor- Fourier modes to l£ and I~, 
respectively. The identity operator is I = 1+ + I~. 

Furthermore, consider the same sequence of positive numbers A n given in (j3.11|) . 

4.4. Analyticity improvement. 
Lemma 4.3. If 5 > 0, 

Pn<^-^ and r n >-^-(p' n + 5)\\ T T^- l l (4.34) 
An-i log 2 

then C n as a map from — E)^4 Pri _ ljr nA 7n to (I — E)A' p , r , is continuous and compact 
with 



\£n\\< (l+ 2 J + 7T^r^) — ' ( 435 ) 

1 5 2r' z \og2J ji n 



Proof. Let F G - ^)A Pn , ur R A 7 „, 

E = {(0, v): u & (N X {0}) d } and J n = {k G Z d : \k ■ w (n) | < <r„||fe||}. (4.36) 
Using Lemma f4. II and (|4.28J) we have 

^ n = fin\\ T T {n) \\ r' + \\b n (H)\\ <- + -\\Qn l -i\\ \WA < ~- (4.37) 

4 r 2 
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We want to find an upper bound on 

\\F oL n {-^ n -)\\'^ >rl 

< (1 + 2^11^-^11 +/i„|| T T(") || \\u\\/r>) iF^l^e^ll^" 1 ^ 

Cl-E 

^ E fl + W - lfc H + ^^e^H-ll) l^l^rHe^H^" 1 ^, (4.38) 
lU-E V 4r 4n / 

where we have used the inequality Ce _<5 ^ < S^ 1 with £ > and again a choice of \x n 
verifying (|4.28|) . Here £ n = | \og(r/ip n ) > |log2. 

Consider separately the two cases corresponding to the definition of the resonance 
cone 

We deal first with the modes corresponding to k G J„_i — {0}. By (j3.11j) and ()4.34j) 
each of these modes in ()4.38j) is bounded from above by 

1 + ^F + ^- ) r lk " ePn_1 " fe "- ( 439 ) 



5 ' 2r /2 log2^ 
Now, consider \\u\\ > T n ||fe|| with k 7^ 0, so that 

|| T T( n ) _1 fe|| < r- 1 !! T T< n ) _1 || (4.40) 
These modes in ()4.38j) are estimated by 

+ 5 + 4r' 2 £ n J V J S \ i + 5 + Ir'Hogl) ' 

(4.41) 

where we have used (|4.34|) . 
Finally, we get 

ll F °^(^OII^< + y + H^llp»-i.r, 

and (jOSj) follows from (jOTJl . 

The compacticity of the operator follows in the same way as in the proof of Lemma 

EH □ 

Let < p" n < p' n and the inclusion 

l n . A' p , n y -> A' p ,, y , H ^ H\V p ,, y . (4.42) 

The norm of the k 7^ modes can be improved by the application of such inclusion. 
That is done in a similiar way as in Lemma 13.61 therefore we omit here the proof and 
just repeat the result for convenience. 

Lemma 4.4. If 4> n > 1 and 

O<p^< P ;-log(0 n ), (4.43) 

then 

||J„(I-E)|| <0„-\ (4.44) 
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4.5. Elimination of far from resonance modes. The theorem below states the 
existence of a symplectomorphism isotopic to the identity that cancels the far from 
resonance modes of a Hamiltonian close to H®. 

For given p n ,r',e,u > 0, denote by V £ the open ball in A' +l/r , centred at H® with 
radius e. 



Theorem 4.5. Let < r < r' , a n > 2r'\\Q n \\ and 



48||#o|| r /(27r + 1) 2 (1 + 2vr 



For every H G V £n there exist an analytic symplectomorphism g: T> Pnr 
Af n r satisfying 



KHog = 0. 



This defines the maps 
and 

which are analytic, and verify 

\mH)-id\\ pntr 



H 
H v- 



A 



2d 



I + A 
Hog, 



(4.45) 

(4.46) 
(4.47) 

(4.48) 



< 



\lZH\ 



Pn,r 



rrUII < - 



\h-h: 



On/ 

n\\p n +u,r' • 



(4.49) 



C 



i>2d 



(4.50) 



\\U{H) 

If H is real- analytic, then <3(H)(R 2d ) 
We prove this theorem in Section \A. 21 

Recall the definition of the translation R z on C 2d given by (jH.25j) . 
Lemma 4.6. In the conditions of Theorem \4.5\ if x G IR d and H G V En , then 

<8(HoR x ) = o<S(H) oR x 
on V Pn , r . 

Proof. If g = (5(H) is a solution of (jQB)) in V Pnir , then g = R x l o®(H)o R x solves the 
same equation for H = H o R X1 i.e. 1~H o g = in T> Pn ^. □ 

4.6. Convergence of renormalization. For a resonance cone and p n > 0, the nth 

step renormalization operator is defined to be 

T^n =U n ol n o C n o lZ n „x and TZ = U Q , 

where lA n is as in Theorem 14. 5[ cancelling all the I~ modes. Notice that if H + (y) = 
H°(y) + v-y with H° non-degenerate, IZ n (H + ) = H® for every v G C d . This means that 
the renormalizations eliminate the direction corresponding to linear terms in y. From 
the previous sections the map IZ n on its domain of validity is analytic by construction. 
In addition, whenever a Hamiltonian H is real-analytic, the same is true for TZ n (H). 



34 K. KHANIN, J. LOPES DIAS, AND J. MARKLOF 

Let r' > r > 0, p > and fix a sequence a n such that er > 2r'||Q||. To complete the 
specification of the resonant modes and of e n in Theorem 14 .5\ take 

according to Lemma FQ1 Consider also the constants v and S as they appear in Theorem 
14.51 and Lemma f4. 31 respectively. For n G N, we define the non-increasing sequence 

a 2 A 2 6 C 2 



9 w = min Q n _ x " J] 



(4r'||Q||) 2 JL = 1 ||TW|| 2 || T rW|| 2 ' 



£n TT ^If'W i Ui, (4.52) 



nr=i ii rw-1 ii 3 'y 2io c 3 n^ (i) ~ 1 ii 2 ii t tw _i i 

with Go = 1. The sequence p n is as defined in (13. 44 J) with 

'l2||#°||r' \ \ Vn\Qn -l 

( G \ 1/2 

Pn 

Here. 




> 1, 



(4.53) 



2 8 C„max{l, \r] n \}Q 




c<= 1+ ^nnr l + _ - ||T r W||3 >1 . 



Notice that n is our choice for Lemma f4. 41 

Define the function i— > <6(u;) as in ()3.46j) but using the present choice of <p n . 

Theorem 4.7. Suppose that det(Q) ^ 0, 

B(u>) < +oo, (4.54) 

and p > B(uj) + v . There exists c,K>0 such that if H G A p y and \\H — H°\\ p y < c, 
then H is in the domain ofTZ n and 

\\n n (H)-n n (H )\\ Pn:r <KQ n \\H-H°\\ p y, nGNU{0}. (4.55) 

Proof. Let £ > and p = p — v — £ >0 such that p > B(uj). Hence, by the definition 
()3.44|) of p n , there is R > satisfying p n > RB~\ for all n G N. 
If c < e we use Theorem 14.51 to get Hq(H) G lQ^4 P0jr with 

||ft (#) -n Q (H°)\\ pihr < KQ \\H - H°\\ py 

for some K > 0. 

Now, suppose that i? n _i = TZ n -i(H) G I^_i^4.p n _i,D ^ G N, and 

H-Hn-l — ^n-lllpn-i,r < ^©n-l||# ~ # lUr'j 

n-1 

iiQn-iii<iiQiin^i^iii TW n ir rW n. 



2' 

i=l 
n-1 



(4.56) 



IIQnMl < IIQ-III 2 ^ 1 !^!" 1 !!^" 1 !! ir r(0_1 l 



1=1 
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So, for c small enough, we get 

rsl/2 n-1 nSA/^Wrpfi)- 1 ^ II Tr-fi)- 1 !! „2 

11^.11 «fen 2C " r |J 2 " r " * 32^- («7) 

Thus, Lemma HH] is valid and as a consequence ||6 n (fT n _i)|| < r/8. 

After performing the operators and Z n , we want to estimate the norm of the result- 
ing Hamiltonians. The constant and non-constant Fourier modes are dealt separately 
in 

l n C n (H) = H° n + C n (EH n „t - Hl_ x )+X n C n {{l - E)H n _i). (4.58) 

For the former we use Lemma 14.21 and for the latter Lemmas 14 .HI and 14.41 That is, the 
definition of \i n implies that fM n < 8r ,^ T T ( n )^ and 

\\C n (EH n ^ - H^t, < 2 7 K (l + -L) f 7 ^) » 2 n \ Vn \ || T T<»)|| i »e fV _ 1 I| J ff - H%s 

K 

< ; ©t) \ \H — H \\or'- 

~ 2(1 + ^l2\\H%,/e n ) 

(4.59) 

Furthermore, (f) n yields 

~ / 2n r \ 

\\l n C n ((I - E)Hn-l)\\' p ^ <K[l + - + 2r ,2 log2 ) ^K^n^n-lWH ~ H°\ 

K 

2(1 + 0lpSV^) 



< ~rr. / „„„ = — ; Qn\\H — H°\\ p>r i. 



(4.60) 

Moreover, assuming c to be small enough, we obtain from ()4.22|) . ||<5n-i|| _1 < HQn-ill 
and (P37jl that 

IIQnll < vn\vn\ \\t^\\ || T TW||||Q n _ 1 ||(i + ler-^e^xiig^ir 1 ) 

<iMin^Nr«iiii T T«ii<^. (4 - 61) 

i=l 

By using (|4.25jl and again (|4.57jl . 

| T (n)-i|| || ttW^iiiiq^j 



H n \ Vn \(l - lQr-icKQ^WQ-Wl) 

n 

^IIQ-III 2 ^ 1 ^!" 1 ^ -1 !! ll Tr(<rl l 



(4.62) 



i=i 

The Hamiltonian I n C n (H n ^i) is inside the domain of U n since for c small enough 
\cKQ n < e n and ||Q n || < a n /(2r'). The result follows from ijOSjl . □ 

Cf. Remark 13.1 II on how to generalise the above for a small analyticity radius p. 

Lemma 4.8. If oj = ( f ) G M. d is diophantine, then (|4.54j) is verified. 

Proof. The proof follows the same lines as for Lemma \3.12\ using the same choices 
of a n and t n . In fact, to show (|4.54jl it is only necessary to check that the series 
£S n |logk+i||, E5nlog||r(" +1 )||, ££ n log||u;(" +1 )|| and £ B n | log a n+1 | converge. 
This is already done in the proof of the lemma. □ 
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4.7. Construction of the invariant torus. In the following we will always assume 
to be in the conditions of section 14.61 We use Theorem 14.71 to determine the existence 
of an o;-invariant torus for the flow of analytic Hamiltonians H close enough to H° 
(Theorem II. 2|) . This follows from the construction of an analytic conjugacy between 
the linear flow on T d of rotation vector u> and an orbit of H . We will use the notations 

n n 

K = Y[vi and Xn = JJ^- 

i=l i=l 

Let the set A be given by 

A = {H G A p , r > :\\H- H°\\ Pt r> < c}, (4.63) 

which is contained in the domain of TZ n for all n G N U {0}. Given H G A, denote by 
H n = TZ n (H) G I^A Pntr . It is simple to check that 

H n = —[(1- E )H o go o o 9l o • ■ ■ o I£] o g n 

Xn 

= —((1- E )H o 9o o \V X {H) o gi o P 1 {H)~ 1 ]o ( 4 - 64 ) 

Xn 

■ • • o [Vn-xiH) o g n ^ o Vn-xiHY 1 ] o V n {H)) o g n . 
Here, g k = <3 k (C k (H k -i)) is given by Theorem 14.51 at the kth step and 

L k : (as, y) i-> L k (x, fx k y), (4.65) 
with L fe as in ()4.8|) obtained for H k -\. In addition, the affine symplectic map V n (H) = 

V n (H): (x,y) i ► (pW" 1 !Bl $i(ff)...$„(ffn-i)(y)) (4-66) 
and P (ff) = Id, where $* (i/) = fi k T T^y + b k {H k _ x ). Notice that 

. . . $ n (# n _i)(?/) = Xn T P {n) y + v n (H), (4.67) 

with 

n 

v n (H) = h(H) + T P {l - 1) k{H l _ 1 ). 

i=2 

Let 



a n (H)= lim $„(#„_!)... $ m (# m _!)(0) 

m— >+oo 



+oo 

To-'Cnl To-ifi— lit / ^ 

i=n+l 

if the limit exists. If that is the case, 



(4.68) 



a(H) = ai (H)= lim u n (i2") (4.69) 

n— >+oo 

and 

a(#) - v n (H) = xn T P {n) a n+1 (H). (4.70) 

Lemma 4.9. For H G A taere zs a n (H) G S r / 2 C C d and the map a n : H a n (H) 
from A into B r / 2 is analytic taking any real-analytic H into M. d . 
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Proof. From Lemma 14.11 we obtain <C 1 for any k G N. Thus, by the 

definition of there is < A < 1 such that 

f ^...^ 1 \\ T TV>... r Tl i -%(H ir . 1 )\\ « A*"", (4.71) 

where 1 < n < i — 1. Hence, f)4.68j) converges and each a n (H) is well-defined in C d , 
unless H is real which gives a n (H) G IR d . The maps if i— > a n (H) are analytic since the 
convergence is uniform. □ 

Lemma 4.10. There is an open ball B centred at H° in A such that we can find 
sequences R n , r n > satisfying R-\ = p, r_i = r' , 

R n + 2tcKQ 2 J 3 \\H - H\ y < P„_i < ^~_\^ , (4.72) 

r n + K&T\\H - H\ y < r n _! < 2|| T *^ r y (473) 

H G B, and 

lim R^&J* = 0. (4.74) 

n— >+oo 

Proof. Let p* = minp n . Since Xn is decreasing, it is enough to check that 

{n n 
Avnii TW ii" i 'X»nn Tr(<) " 1 ii" 1 
i=i t=i 

for some < A < 1 by taking R n = Ci\~ n< df/ 3 and r n = c^On^ 3 with small constants 
C\,C2 > 0. Thus, the inequalities ()4.72|) and f)4.73j) hold whenever we take a sufficiently 
small bound on \\H — H°\\ P)r . The limit ()4.74j) is now immediate. □ 

For each H G A, consider the isotopic to the identity analytic symplectomorphism 

W n (H) = V n (H) o (d^CniH^)) o V^H)- 1 (4.75) 

on V n (H)V Pn r . In particular, W n (H°) = Id. Notice that for H real-analytic, 

W n {H){R 2d ) c R 2d . 

For a given if G A, define the norm ||X|| n = \\X o V a (H)\\R n , r „, whenever X o V a (H) ^ 
A 2 ™ _ , where we have introduced the vertical translation 

(x,y) ^ (x,y + z), (4.76) 

for any z £ C d . 

Lemma 4.11. W n is analytic on B, and satisfies 

W n (H): V a(H) (V pn , rn ) - K( fl )Pp„^„J 

and 

\\W n {H) - Id ||„ < K'Q 2 J 3 \\H - ff°|| Pir , H E B, (4.77) 
/or some constant K' > 0. 

Proo/. For if G A and (cc, y) G Di^, 

||ImP (n) :r|| < ||P (n) \\R n /2n < p n /2n, 

WK'iH^) . . . *?{H)(y + a(H))\\ = fc 1 T P {nr \y + a(H) - v n (H))\\ (4.78) 

< Xn'll T P W ~ 1 ||r r ,+ ||a„ +1 (i?)|| <r. 
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Therefore, V n {H)- x o V a(H) (V Rn ^ n ) C P Pn , r . Moreover, using fOSJ, 

\\W n (H) - Id ||„ = \\V n (H) o [(& n (X n jC n (H n ^)) - Id] o Vn(H)- 1 o V a{H) \\ Rn , rn 

< e^WPnfflW \\l n C n (H n ^) - H% n y (4.79) 

<^Q 2 J 3 ,\\H-H\, r , 

where V n (H) corresponds to the linear part (x,y) \— > (p( n ) 1 x, \ n T P ijl ">y) of V n (H) 

which has norm bounded by \\V n (H)\\ < ||P(") \\ + Xn \\ T P {n) \\- 
Now, for (a;, y) G V Rn)Tn and H e B, 

\\irilmW n (H) o V a{H) (x,y)\\ < \\ ]m(mW n (H) o V a(H) (x,y) - x)\\ + || Imcc|| 

< \\W n {H) - Id || n + R n /2n < R n ^/2n, 
\\n 2 W n (H) o V a{H) (x, y) - a{H)\\ < \\n 2 W n (H) o V a{H) (x, y) - y - a(H)\\ + ||y|| 

< \\W n (H) - Id || n + r n < r„_i. 

So, W n (H) : V a{H) (V Rnir J - ^(2?^,^). □ 
Consider the analytic map T n on 5 satisfying T n (if) : V a{H) (V Rntrn ) -> V a{H) (V p y), 

T n (H) = W (H)o...oW n (H). (4.80) 

We then rewrite ()4.64j) as 

H o T n (H) = ^H n o Vn(H)- 1 + E(H), (4.81) 

where E(H) represents a constant (irrelevant) term. Since each W n (H) is symplectic, 
thus T n (H) is symplectic and H o T n (H) is canonically equivalent to the Hamiltonian 
i? n . In particular, if H n = H® for some n, there is an o;-invariant torus in the phase 
space of H n . We are interested in the general case, H n — H° — > as n — > +oo. 
By a simple adaptation of Lemma lH.15| we have the following result for T n . 

Lemma 4.12. There is c > such that for H £ B 

\\r n (H) - r n ^(H)\\ n < ce 2 J 3 \\H - H°\\ py . 

Consider the Banach space Cp er {R d , C 2d ) of C 1 functions Z d -periodic, endowed with 
the norm 

\\f\\ c i =maxmax|| J D A; /NII- 

Lemma 4.13. There exist C > 0, an open ball B' C B centred at H° and an analytic 
mapT: B' -> Diff per (M d , C 2d ) such that, for every H G B' , T(H) = lim n ^ +00 T n (H)(-, a(H)) 
and 

\\T{H)-{ld,a{H))\\ C i <C\\H-H°\\ py . (4.82) 
IfHeB' is real- analytic, then T(H) G Diff per (M d , R 2d ). 

Proof. For each H G B, by the first inequality in (|4.5jl . 

|[r n ( J ff)-r n _ 1 (i?)](.,a( J ff))|| c71 < max sup ||D fc [T n (ff)(aj, a(F)) - T n ^(H)(x, a(H))]\\ 



k < 1 xeD n 



tin 

(4.83) 

which is estimated using (|4.74jl . Hence, T n (H)(-,a(H)) converges in the Banach space 
Cp er (R d , C 2d ), and ()4.82|) holds. The convergence of T n is uniform in B, thus T is 
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analytic. If H is sufficiently close to H°, T(H) is in fact a diffeomorphism as the space 
of close to identity diffeomorphisms is closed for the C 1 norm. Finally, for H real- 
analytic we have T(H)(M. d ) C M. 2d in view of the similar property for each W n (H). □ 

The Hamiltonian vector field of H is denoted by Xh = JV-ff, where J: (x,y) i— > 

(y, ~ x )- 

Lemma 4.14. For H e B' , we have on M. d 

X H oT(H) = D(T(H))u>. (4.84) 
Proof. Since T n (H) is a symplectomorphism, we have for x £ M d , 

y n (aj) = o T n (H) o 0) - D(V n (H)) o 0) Xflo(aj, 0) 

= [D(V n (H)) o K(^)^or„(^)oy a(H) -i/o](a;,0). 
Hence, 

||y n (aj)|| < \\D(T n (H))(x,a(H))\\ \\V[H o T n (H) o V a {H) — H°}(x, 0)||. (4.86) 
We can now estimate the above norm by recalling (|4.81|) in 

V[H o T n (H) o V a (H) - H°] =^V[(H n - H° n ) o V n {H)~ l o V a(H) 

+ El o Vn(H)- 1 o V a{H) - Hi o ^(tfor 1 ] 

and 



(4.85) 



(4.87) 



^V[H° n o P n (tf )-* o V a(H) - Hi o V n (H yi] 

n-l 

T| 



[a(F) - v n (H)] 
Notice that by induction we get 



Q + J2 ^-pwV^^)) TP®" 1 



(4.88) 



n-l 



(4.91) 



Qn = XnXnP H Q T P<"> + T P®'* & F® {^{Hi)) TP®"* T P W . (4.89) 

since ESfolAiD-lpW" 1 !! II T P (?rl ||e ? , « 1, 

l|^||V[(F n - P°) o V n (Hy l o V^hjK*, 0)|| « V^ 1 !! ||F n - P°|| Pn , r (4.90) 
and 

||a(P) - <P)|| < Xn|| T P (n) || IK+l(P)|| 

«Xnll T P (n) l|e ri 

we have that 

\\V[H o F n (H) o - H°](x, 0)|| « e^ 3 . (4.92) 
Finally, from the convergence of T n and 

\\DT n (H)(x,a{H))\\ « -^-\\T n (H)\\ n « -L (4.93) 

we find that ||Y^(a;)|| converges uniformly to as n — > +00. □ 
Lemma 4.15. If H E B' and x E R d , then 

T(H o i^) = P^ 1 o T(P) o ^ (4.94) 
where R x : z ^ z + x is a translation on C d . 
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Proof. For each n G N, (|4.32|) implies that V n (H o R z ) = V n (H) and we know that 
V n {H) o R-\ n)z = R- 1 o V n (H), z G C d . So, from LemmaiH 

w n (H o r x ) = v n (H) o e n (c n n n ^(H o r x )) o Vn(H)- 1 

(4.95) 

= R- 1 oW n (H)oR x . 

Thus, T n (HoR x ) = R' 1 oT n (H ) oi?^ and ()4.94|) follows from the convergence of T n . □ 
The flow generated by X H is denoted by <$ H taken at time t > 0. Hence, 



if 



We prove below the existence of an invariant torus T for H, i.e. an analytic conjugacy 
between </>#|t and i?^. 

Theorem 4.16. Let D C R d be an open ball about the origin. If H G C"(T d x £)) 
is sufficiently close to H°, then there exist a G M. d and a C w -diffeomorphism 7: T d — ► 

4°7 = 7° flwt, < > 0, (4.96) 
and T = 7(T d ) ~ T d zs a submanifold homotopic to {y = a}. Furthermore, the maps 
H 1— > a and i7 1— > 7 are analytic. 

Proof. The lift to M d x £) of is assumed to have a unique analytic extension to T> p y. 
Consider the real-analytic Hamiltonian G = H G A p y. Suppose that G is close enough 
to H° such that G G B' and G o R z e B' for 77 > and z G D n . Then, a = a(G) and 
7 = r(G)|[ 01 )d, C 1 and homotopic to (Id, a), verifies ()4.96|) . This follows from (|4.84|) 
and the equivalent equation 



d 
It 



d 



(^°7) = t (7°U 
f=o dt f=o 

which we integrate for initial condition 4>° H = Rq = Id. 

We now want to extend analytically 7 to a complex neighbourhood of its domain. 
Take 7(2) = R z o T(G o R z )(0), z G D v . The maps z ^ G o R z and H i-» r(#) 
are analytic and Cp er (lR d , C 2d ) 9 g 1— > a(0) is bounded. As 7: — > C M involves their 
composition, it is analytic and Z d -periodic. From ([4.94)1 . for any x G R , we have 

7(aj) = r(G) o ^(0) = r(G?)(o:) = 7 (aj). 
Finally, since T is analytic, the same is true for the map H 1— ► 7. □ 

Remark 4.17. A quasiperiodic invariant torus T is always Lagrangian. If in addition 
T is homotopic to {y = a}, it is the graph of a function ip G C UJ (T d , D) (see [T2]). 

Appendix A. Elimination of modes 

A.l. Homotopy method for vector fields. In this section we prove Theorem 13 . 71 
using the homotopy method (cf. [20] )• As n is fixed, we will drop it from our notations. 
In addition we write p' = p n and p = p n + v. We will be using the symbol D x for the 
derivative with respect to x. 

Firstly, we include a technical lemma that will be used in the following. 

Lemma A.l. Let f G A'. IfU = ld+(u, 0) where u: D p i x B r — > D( p _ p iy 2 is in A p ' jr 
and ||n|| p ',r < (p — p')/^ , then 

• ||/<>l7||p#, r < 11/11 {p+p')/2,r, 

. \\D x foU\\<\\f\\[ p+pl)/ ^, 
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• 11/ O U - /|| p / ir < \\f\\{ p+p/)/2 ,r \\ U \\ P ',r, 

. \\D x f o U- D x f\\ < -^\\f\\' p , r ||u||,, r . 

The proof of these inequalities is straightforward and thus will be omitted. Now, 
assume that 

5 = 42e/a < 1/2. 

For vector fields in the form X = uj + 7r 2 + /, where tt 2 : (x, y) i— > y is seen as a function 
in A', consider / to be in the open ball in A' centred at the origin with radius e. 
The coordinate transformation U is written as U = Id+(w, 0), with u in 

B = [u e I~A' p , tr : u: Dp, x B r ||w||p/ ir < 5} . 

Notice that we have 

ru*(x) = r{Duy\u + 7T 2 + / oc/,o) 

= (r(J + J D a! «)- 1 (a; + 7r 2 + /oC/) ) 0). 

From now on the parameter r is omitted whenever there is no ambiguity. Define the 
operator F: B — > 

F{u) = r(I + D x u)-\u + ir 2 + foU). (A.l) 

F(u) takes real values for real arguments whenever u has that property. It is easy to 
see that the derivative of F at u is the linear map from I - -4.|y to l~A p i: 

DF(u)h = r(I + D x u)- 1 [D x f oUh (A.2) 
-D x h (I + D x u)- 1 (w + tt 2 + / o U)\. 

We want to find a solution of 

F(u t ) = (l-t)F(u ), (A.3) 

with < t < 1 and "initial" condition no = 0. Differentiating the above equation with 
respect to t, we get 

DF(u t )^ = -F(0). (A.4) 

Proposition A.2. If u G i3, i/ien DF(u)" 1 is a bounded linear operator from I~A P > to 
I~A' p , and 

WDFiu)'^ < 5/e. 

From the above proposition (to be proved in Section IA.1.1|) we integrate (jA.4|) with 
respect to t, obtaining the integral equation: 



/ DF(u 8 )- 1 F(0)da. 
Jo 



u t = - DF{u s )- l F{Q)ds. (A.5) 



In order to check that Ut G B for any < t < 1, we estimate its norm: 



\u t \\'< tsuppF^)- 1 ^)!!' 



< tsup||L>F( U )- 1 ||||I-/|| / ,<M||/|| l0 V£, 

so, ||wt||'/ < 5. Therefore, the solution of (jA.3|) exists in B and is given by (jA.5|) 
Moreover, if X is real-analytic, then ut takes real values for real arguments. 
It is now easy to see that 

U*X - X° = 1+ ^2(-D(U t - ld)) n X° + l + U* t f + (1 - t)rf. 

n>2 
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So, using Lemma fA. 11 



\\u;x-x%, < 



1 



p' 2 + II./ II,-' 



< 



< 



1 f5 2 \\oj 



(i-t)\\fV 

P '/e 2 + l) \\f\\ P + (1-011/11 



1-6 \ e 



1 + 1-t 



Moreover, \\U;X -X°-I+ f -{l-t)I- f\\ p , = 0(\\f\\ 2 p ), hence the derivative of X h-> C/*X 
atX°isI-tr. 

A. 1 . 1 . Proof of Proposition \A.£H 
Lemma A. 3. If \\f\\'<e< a/6, then 

is continuous and 

i 



\\DF(0)- l \\ < 



*-nf\\ P 

Proof. From ()A.2j) one has 

DF(0)h =r(f-D„)h 

'l-I-fD^Dvh, 



where / h = Df h — Dh f and D w h = D x h (u> + 7r 2 ). Thus, the inverse of this operator, 
if it exists, is given by 



The inverse of is the linear map from I A p > to I A ; a ,: 

D„ l g(x,y) = 



v -p 1 

9k(y) 



2-Kikx 



feel- 



2mk ■ X°(y) 



and is well-defined since Lemma 13.21 implies that 

\k.X°(y)\>a\\k\\/2, 

with k e I~ and y £ B r . So, 

! \^ l + 27T||fc|| 

A* P Ly = >. o sup 



9k{y) 



k ■ X«{y) 



J\\k\\ 



< ^ l±MM|| fft ||^iH 

feez- 
3 

< - 



7TO" 



Hence, H-D^/H < 3/cr. It is possible to bound from above the norm of /: I A' p , — ► ^4 P ' 
by ||7|| <2||/||;,. Therefore, 

ir/^n < -ii/ii^ < i, 
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and 



-i 



-i 



< 



a 



^-6||/L 



The statement of the lemma is now immediate. 



□ 



Lemma A. 4. Given u G B, the linear operator DF(u) — DF(0) mapping I A' p , into 
I~A P >, is bounded and 



u 



4tt 4-2|H[ ^ 

i ~ ' -i ii ii/ 



■u 



\\DF(u) - DF(Q)\\ < 

_ P-P - ii-np' / ■ 

Proof. The formula (jA.2j) gives 

[DF(u) - DF(0)] h = l-(I + D x u)- l [D x foUh-(I + D x u)D x fh 

-D x h (I + D x u)-\u + n 2 + f)oU 
+ (I + D x u)D x h (u> + 7T 2 + /)] 
= T(I + D x u)- l {A + B + C}, 

where 

A = [D x foU-D x f-D x uD x f] h 
B = D x u D x h (u + tt 2 + /) 

C = -D x h (I + Dg-u)" 1 [foU-f-D x u(uj + n 2 + /)] . 
Using Lemma [A. 11 



+ 7T 2 



< 



< (lk+7T 2 | 



47T 



< 



1 



pIMU + 



p'IMIp' 



p'J iFlIp' 



■(/, 



r 01/11 W)/2lHU + IHrp'(ll w + ^||r- + ||/||pO] 



!->'■ 



To conclude the proof of Proposition IA.21 notice that 
\\DF{u)- x \\ < {WDFiOyT'-WDF^-DFm)- 1 



□ 



< 



o 5 



47T 



P~ P 



4-25 
1-5 



e + 



2-5 



OJ + 7T 2 



< 



The last inequality is true if 



e < 5 



a 



25 



U> + 7T 2 



1 + 25 + 



4tt 4-25 
+ 



\-5\p-pl 1-5 



with a positive numerator A" and denominator D in the r.h.s. This is so for our choices 
of e and 5 < |, by observing that 



Hlr = sup ||y|| < a ^^I^IIM^H + 6||w|| < -||w| 
ye B r 2 



(A.6) 
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thus 1 1 + 7r 2 || r < |||w|| and 

25 „ , M1 , a 

U) + 7T 2 || r < 125||o>|| < 



(1-5)2" 11 6 ' 

So, iV > 5tr/6, D <7, and finally £ < < § < N/£>. 

A. 2. Elimination of modes for Hamiltonians. Here we present a proof of Theorem 
14.51 It is similar to related methods appearing in e.g. ^21 Q]. As we have fixed n, we 
will not include it in our notations. 

Let R = R 2 ) and R' = (R[, R' 2 ) be such that R> R' > componentwise. We will 
be interested on the set Q R > of analytic symplectomorphisms g: T> R > — ► T> R satisfying 
g-lde A R , and 

\\g - Id \\ R > <5 = min{(i2i - R' 1 )/2tt,R 2 - R' 2 }. 

We use the notation {•, •} for the usual Poisson bracket associated to J: (x,y) i— > 
(y,—x). In the following R — 5 stands for i? — 5(1,1) and 7r2 : {x,y) i — > y is the 
projection on the second component. 

Lemma A. 5. let < £ < ~. If G E A' R , and \\G\\' RI < £5/(2n + 1), then there is a 
unique analytic symplectomorphism g : T>r>-2S —> < C 2d such that \\g — Id ||jj/_2<5 < £<5 awei 

g = ld+SVGog J (A. 7) 

where g(x, y) = (x, ir 2 g(x, y)), (x, y) G Pr/_25- Moreover, for any H e Ar> 

\\H o g\\ R >-2S < \\H\\r> 
\\Hog-H\\ R ,„ 2S <2Z\\H\\ R , (A.8) 
\\Hog-H- {H, G}\\ R ,_ 25 < 2e\\H\\ RI 
and the maps G i— ► g and G *—>■ H o g are analytic. 

Proof. Define the map T : g i— > Id +JVG o g on the open ball .B in G R >-2S centred at the 
identity and with radius £5. It is simple to check that T(B) C B, in particular T(g) for 
g G -B is symplectic. We now show that T is a contraction on B and thus its unique 
fixed point is the map we are looking for. In fact, whenever g e B we obtain 

\\DT(g)\\ < \\DVGog\\ R ,_ 2S < \\DVG\\ R ,_ S 

< ^iivgii, < ^iicir, < ( . (A ' 9) 

For the estimates in (jA.8|) (the first is now immediate) we introduce the differentiable 
function 

/: {^C: \z\<(}^A R , 

z^Ho(ld +zIVG(Id +z(g - Id))) 

where ( = l/£ > 2. Cauchy's integral formula yields that 

\\Hog-H\\ Rf - 2S = ||/(1) -/(0)||j*_ 2 , 

" 2tt Xl=C 1^-1)1 (A.H) 

<-i-SU P ||/(z)||^- 25 <2£||if||tf. 

C, — 1 |z|=£ 



(A.10) 
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and 



\Hog-H-{H, G}\\ R ^ 2S = ||/(1) - /(0) - /'(0) \\ R ,-26 



2nJ lzl=( |z 2 (z-l)| (A.12) 

(,{(,- l) |z|=£ 

By the implicit function theorem the maps G i— > g and G h- > H o g are analytic. □ 
Lemma A. 6. Let a > 2R2WQW, e' > and H e .4.^ suc/i i/icrf 

||/r-ff°||fl <£' < - ^— — — . (A.13) 

11 UR ~ (27r + l)[l + 2 7 r + (r + l)/ J R 2 ] V ; 

Then there is G G I~*<4/ R , such that 

r(H + {H,G}) = and \\G\\' R , < —J— \\rH\\ R ,. (A.14) 

Moreover, the map H G is analytic. 

Proof. Consider the linear operator associated to H: 

T: TA' R , -> I' An, K i-> I~ (A.15) 

It is well-defined since 

\\r{H,K}\\ R , < \\VH\\ R r\\VK\\ R , 
We will show that : I~A R i — > I~*4/#/ is bounded and 

We start by decomposing any Hamiltonian H = H° + F as 

H{x,y) = Y,H k {y)e^ ikx with # fe (y) = Y,H k , v y". 

k v 

Write D = V2H • Vi, with Vi and V2 standing for the derivatives with respect to x 
and y. The definition of T in (|A.15|) yields 

T{K) = r (F - D ) K = - (l - rFD ^ D K, 

where F(K) = {F, K}. If the inverse of T exists is given by 

T~ l = -D G l (l - rFDo 1 ) 1 . (A.17) 

The map Dq 1 : I~A R > — > I~A' R , is linear and given by 

kez d -{o} v yy>l 
For each ke I~, using (jj^Hj) and < a/(2R 2 ) thus |fc • Qy/k ■ u\ < 1/2, 

W k {y) W k (y) ^ / fc-Qr/ \ n ^ (A lg) 



fc-w(l + ^K) fe w fe-W 
V fc-tt) / n>0 x 
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we get the estimate 



k ■ V 2 #° 



^EE 

R 2 n>0 v 



a 



n+l I 



k 



< 



EE 

n>0 v 
2 



\W k ,„\Bl l 



a 



k\ 



cr\\k\ 
W k \\ R2 . 



(A.19) 



Similarly, we find the bound 
V 2 W k 



k ■ WoH° 



< 



H 2 



2||i/|| \W k , v \Bl 
a\\k\ 



h>\\-i 



< 



2r 

~^R->. 



\W k \ 



R 2 - 



Finally, 



W k Qk 



[k ■ V 2 H ) 2 
It is now immediate to see that 



< 



r 2 <rRa\\k\ 



\Wi 



k\\R 2 - 



(A.20) 



(A.21) 



\Dv l W\\ R < ^-WWU, and HViGDo^WOlU < -\\W\\b>. 

Z7TCT" a 



Moreover, 



W k (y) 



V 2 W k {y) 



W k {y)Qk 



which implies 



Hence, 



k-\7 2 H°(y)J k-V 2 H°(y) (k-V 2 H°(y)) 2 
\\V 2 (D^W)\\ R/ < I+±\\W\\ R ,. 

7YO~K 2 
I T + 1 



(A.22) 



|L> - 1 ||<-(1 + — 
a 



2tt 2nR^ 

As /•': J~A' a -> with ||F|| < 2 ||VF|| fl , < 2^±I||F|| fi (by Cauchy's estimate), 



a 



+ 2vr + 27T.R2 



iVFUfl/ < 1 



and 



-1 







< 









+ 2vr + 27T.R2 



IVFI 



R! 



Thus exists given by (|A.17jl and the estimate (|A.16J) on its norm follows immedi- 
ately. 

A solution of (|A.14|) is simply given by G = I~H). Therefore, < 

H^- 1 !! ... ^ 

Consider the pairs R = (p n + u,r') and R' = (p n ,r), a > 2r'\\Q\\ and H = H as 
given in Theorem 14.51 We are going to iterate the procedure indicated in the previous 
lemmas. Let a sequence of Hamiltonians be given by 

H k = H k ^ o g k , keN, 

where G k and g k are determined for H k _ 1 by Lemmas IA.6I and IA.51 respectively. In 
addition, denote by 

fto-oj* (A. 23) 



,(*) 
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the composition of all symplectomorphisms up to the fcth-step so that Hk — Ho g( k \ 
In order to determine the right domains of H k , Gk and gk, define the sequences 



i=l 



R k = R k ^-A5 k = R-AJ2 S 
R'k — Rk-i — $k, 

with R = R and 



From now on, assume that 



Lemma A. 7. If for every fceN, ||I Hk-i^R^ < e'/2 and 

(2tt + l)e' 
i/ien Qk{T>R k ) C X>r a _ 1 and 

ii^-idik^^lp-^ik^ 

i=i 

Proof. Recall Lemma fA. 51 for £ = ||I~iiffc_i||^ fc _ :L /£ / and check that 

Ibfe - Id |k < \\g k - H lk-24 < -^|l r #fc-ilk-i 
and Rk + Sk < Rk-i componentwise. Now, 

k k 



l " ) - Id = E 9i ° • • ■ ° 9k ~ E 9i ° ' ' ' ° Sk ~ Id 



i=l i=2 
fc-1 



8=1 



Thus, 



(A.24) 



4= 2^ min { 1 '27' r '- r }-^ (A ' 25) 



^ndn^ ^ ^ , (A.26) 



(A.27) 



E(* °---°9k- 9i+l o ■ ■ ■ o g k ) + g k -ld (A.28) 

1=1 
fc-1 

E(* ~~ Id ) ° ° ' ' ' ° 9k + gk - Id ■ 



\\9 {k) ~ M Ik < E " Id II* ^ E ^P^-ilk-!- ( A - 29 ) 
j=i i=i 

Furthermore, as 

g (k) _ g (k-l) = {g (k-i) _ w) Q gk _ {g (k-i) _ M) + {gk _ M) {AM) 
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we get 

\\9 {k) - 9 ik - 1} \\R k <(\\Dg^ - J|k + 1) \\9k - Id II a 



8h r „ / 27T + 1 > 



^p-^lk^l^^plir^ll^ + i, (A . 31) 



□ 



Notice that since e' < ^||if°|k we have 



2 

£' < - e' < \\H\\ R < \\H°\\ R + e' (A.32) 

and also 

~\\H \\ R <\\H\\ R <~\\H°\\ R . (A.33) 
Lemma A.8. For any k eN, if \\I~H\\ R < e' 2 /(8\\H\\ R ), then 

\\l-H k \\ Rk < {^) 2k ~\\TH\\t< £ -, (A.34) 

\\H h - H k ^\\ Rk < ^^lir^-ill^, (A.35) 
||#klk < 2||F|| H . (A.36) 

Proof. We will prove the above inequalities by induction. The generating Hamiltonian 
Gi given by Lemma fA.6l and the symplectomorphism gi by Lemma fA .51 satisfy ||Gi|k < 
5 1 ||I-^||/[(27r + ly], \\ 9l - Id ||^_ Ml < Wl-HW^/e* and l~H x = I~H o ^ - r(# + 
Hence, 

p-#i|k < \\H o 9l -H k - {H, Gx}\\ Rl < 2 r 111 ||#|k (A.37) 

and 

Htfr-tflk < HV^Ikll^i-Wlk < ^^Ill-^II^HiflU < pUFfflUllfflU. (A.38) 

Thus, |A~34jl and fA~35j) are valid for k = 1 and so is ()A.36jl because ||i?i|k — ll-^i ~~ 
H\\ Rl + \\H\\ R . 

Now, assume that the inequalities are true for k. Under these conditions, Lemma [A. 61 
guarantees the existence of Gk+i so that 

HGWik +1 < j2^W' U ' HkllRk (A ' 39) 
and Lemma IA.5I yields g k +i- Therefore, \~H k+ \ = l~H k o g k+1 — I~(Hk + {Hk, Gk+i}) 
and 

P~#fc+i|k +1 < \\H k — H k — {H k , Gk+i}\\n k+1 

21 ? ) » H4r * (A.40) 
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Similarly, 

\\H k+1 - H k \\ Rk+1 <\\VH k \\ Rk+1 \\g k+1 - Id \\ Rk+1 



4 . 
< — ||I H k \\ Rk \\H\\ R . 



Finally, making use of the above inequality, 

fc+i 

H#jfc+ilU fc+ i < \\H\\r + 2J - •Hi-illflfc+j 



i=i 

fc+i 



< + -^£l|rir«-i|k_ 

i=i 



< \\H\\ R +\\H\\ R J2{ 

i— 1 \ 



' i //|,,|: //|/,. x 2 1 



(A.42) 



i=l 

k 



e' 2 



8=1 



□ 



Theorem 14.51 will now be a consequence of the following result. 
Theorem A. 9. If 

\\H-H°\\ R <e= „ - < „ „ , (A.43) 
11 llR 12\\H°\\ R ~ 8\\H\\ R ' V } 

then there exists g = linifc^ +00 E G R > such that T~H o g = on T> R i. Furthermore, 
the maps <5: H 1— ► g and U: H 1— > H o g are analytic, and 

lb -W Ik < -p-ffHfl (A.44) 



|^o^-iJ || R , < |i + A / 12 H g °llM \\H-H°\\ R . (A.45) 

Proof. Lemmas IA.7I and I A. 81 imply that the sequence converges to a map g: T> R i — > 
T> R which is analytic and symplectic, and = lim k _ >+OQ H k = H o g. Moreover, 
I~Hog = I^Hoo = 0. Since the convergence is uniform, the maps H 1— ► g and H \— > Hog 
are analytic. 
Notice that 

~" 'AWHWrWI-HWrY 1 ' 1 4||#|U|lI-#|| iJ ^ /4||F|| jR ||I- J ff|| x 



El ^W 12 WRW 11 n \\R \ < ^\\ n \\R\\^ n \\R , \ ^ 
1 e' 2 J ~ e' 2 ^ V 

i=l x ' i=l 



< {l + l^B-H\\ R ) 4 A\I-H\\ R (A.46) 



3£ ,2 ir -nicy £/2 

<20jjff|U 1 

- 3£/ 2 11 UK - e ii iiit 
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The inequality in (jA.44|) follows by taking the limit k — > +00 in (|A.27|) . That is, 

+00 „ 1 

\\g - id \\ R > < -4ll r ^-ilk-! < -\\ rH \W- (A.47) 



e- e 



Now, 



\Hoo — H°\\ Rk < \\H — H°\\ R + ^ \\Hi — 



R, 



i=l 



< 12^ n 0| 



e 

where we have used Lemma IA.8I and the fact that ||I _ i7||ft<||if — H°\\r. □ 
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